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ABSTRACT 


The purpose of this thesis is the investigation of 
mathematical methods for treating collisions of light nuclei. 

Our understanding of nuclear fonts 1s derived largely 
from the study of simple nuclear systems. Bound nuclei 
consisting of two to four particles have been theoretically 
treated with considerable accuracy and yielded many valuable 
results. Another important source of information are nuclear 
collisions. However, up to the present, only two particle 
collisions have been theoretically treated in a systematic 
and reliable manner. 

This thesis is concerned with collisions involving any 
number of nucleons. The problem is to calculate the scattering 
properties of nuclei from assumed nuclear interactions. These 
theoretical results can then be checked against experimental 
data,thus providing a test for the validity of the interactions 
which were postulated. As the wave function of a nuclear 
n- body system extends over Sn dimensions and is highly 
irregular, on account of the great strength of nuclear forces, 
one cannot hope to determine it with great exactitude for 


systems involving more than two particles. For this reason 


et 


we have developed variational methods which yield the 
scattering properties of nuclei with considerable accuracy 
without requiring an exact knowledge of the wave function. 

The wave function is assumed to satisfy a Schroedinger 
equation in which nuclear forces are described by two-particle, 
velocity-independent potentialse When the colliding nuclei 
are far apart, the form of the wave function is known: It 
consists of a product oe Punetlens describing the internal 
states of the nuclel, and a free particle function, R(r WP he 
the relative position vector e of their centers of gravity. 
To predict the experimental scattering properties of the 
nuclei the detailed characteristics of F(¥) must be 

determined from the assumed interactions. 
If one considers collisions in which the nuclei have 
a definite angular momentum of relative motion, & » the 


-> 
function F(*) has the asymptotic form 
F(*) aay + (sin an + Korn ¢ . Cor An) Se, 


where k is the de Broglie wave number of relative motion, 

) a | ' Sy is the angular function and NL the phase 

shift corresponding to the angular momentum en The phase 

shift Ne determines the scattering cross section. 
Similarly, the function F(P ) corresponding to a 

plane wave of one nucleus incident on another has the 

asymptotic form 


-> aan 
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Here Lk, | ~jicle i, k, is in the direction of incidence 
> > 

and f(k,, k) is the scattering amplitude in the direction 


of k. 

We seek to determine the phase shifts and scattering 
amplitudes, given the interaction between the nucleons. In 
collisions of two nucleons the wave equation separates and 
mts. solution aes eee no essential difficulty. But when the 
colliding nuclei are complex, the Schroedinger equation can 
no longer be separated. Hence an exact solution is not 
possible and one must take recourse to approximation methods. 

We have followed two lines of attack, both of which 
have been successful in bound state problems: The Raylelens 
See variational method of undetermined coefficients and the 


variation-iteration method. 


Rayleigh-Ritz Not henhie 
| The principle of one of the procedures adopted may be 
Llillustrated by the one dimensional Schroedinger equation, 
2 . 
= + Ao - v(x)} u(x) = 
which is to be solved subject to the boundary condition that 
u(O) = 0. The interaction potential V(x) may be supposed to 


vanish for x Za, so that 


\V 
S 


u (x) : Sim Ax + tam y-cor hx ). SC 


We now form the expression 
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° A r | 
T= a(x) | sas +h - Ux) u(x) de 


C 
and consider its first varistion, restric ing cule eee 


functions, u by the conditions uy, (0) = O and 


ey 
Uy (x) — > nme hoe + Kann, wookx , Q> X Doo 


After an integration by parts, we find 


Since for the correct solution u, I =O, it follows that 


tan a determined from the equation 


kK kan » = T +{town 


is stationary, relative to any admissible variation of u. 
When this equation has been re-written in a homogeneous 


morm, one can use a trial function of the type 


on the right hand side and determine the c's from the 
stationary property of a « As in the case of bound 
state problems this leads to a determinantal compatibility 
equation. ‘In contrast to secular equations this is a linear 
equation for the determination of tan " » regardless of the 
number of parameters. 

This and similar methods are applied to the following 


cases; 





A. Two Particle Collisions. 

ls Phase Shifibsy SRiatesso reas expressions are developed 
for the phase shifts in the case of central forces and for 
the proper phase shifts of the scattering matrix in th case 
of tensor forces. Some numerical calculations show that an 
accuracy of a fraction of one per cent is attainable without 
undue labour. 

2. ocagtering Amplitudes. The variational method 
developed permits .an application of the Rayleigh-Ritz method; 
it also leads naturally to the Born approximations and 
eoririiagee'™S variation-iteration formulation ( J.S. Schwinger, 
unpublished lectures ). 

o. Level Width. A stationary expression for the width 


of a resonance level is derived. 


Be Collisions of Composite Nuclei. 

1. Neutron Deuteron Scattering. This case is treated 
in some detail, as a typical example of composite collisions. 
A preliminary calculation is carried out for the scattering 
cross section at zero energy and reasonable agreement 
( within about ten per cent ) is obtained with the experimental 
value of 3.2 barns. | | 

2. General Composite Collisions. A few features of 
the Rayleigh-Ritz method in rather general collisions, which 
may involve excitation or transmutation, are discussed. 
Variational principles for the elements of the scattering 


matrix are derived. 
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Variation Iteration Method. 

Instead of employing trial functions with undetermined 
coefficients in the stationary expressions for the phase shifts 
and scattering amplitudes one can also use iterated functions. 

Thus in the one dimensional problem the n-th iterate 


is defined in terms of the (n-l1)-st by the equation 
oD 


i. 6 (x, x') V (x') U 
where : 


cAy 
Glxx') 2b sin hoe ra hays Sh sinda sinds, 


Professor Schwinger ( J. S. Schwinger, unpublished lectures ) 
has applied this method very successfully in collisions of 
two particles. Since the iteration improves the wave function 
very rapidly this method converges extremely well. Schwinger's 
method is generalized from two particle collisions to 
composite collisions. The persistence of the internal 
interactions in the asymptotic part of the configuration 
space causes considerable difficulty in the construction of 
a suitable Green's function. Of several methods which are 
investigated only one ( involving a Green's function in three 
dimensions ), is found at all feasible. However even this 
procedure is so complicated, that unless a more suitable 
version of the variation-iteration method can be found, the 
Rayleigh-Ritz method appears to be more useful for the treatment 
of composite nuclear collisions. . 7 oY 

ur, be fukngy 

N. Abdote: 
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INTRODUCTION. 


A great deal of information about the nature of nuclear 
forces has been derived from a comparison between experimental 
and theoretical studies of Sunanemeuc eee systems. Among the 
nost important data so compared are the binding energy and the 
magnetic dipole and quadrupole moments of the is nell the 
binding energies of se. ae and He" a: and the scattering 
cross sections of neutrons and protons by auuonsaale 

However, while the binding energies of nuclear S- and 4- 
body systems have been the ov eeieniNl ame ompuced and have led to 
some very important conclusions about nuclear forces (eege the 
forces operative between neutrons ), no very satiecseeers 
scneme for mathematically treating collisions which involve 
more than two elementary particles has so far veen givene In 
what follows some promising methods to solve such problems 
will be investigated and preliminary results presentede 

Theoretical. calculations of collisions of light nuciel 
are worthy of interest for the following reasons: 


(1) They may confirm or weaken our confidence in existing 


nuclear theories; in particular, they may throw light on the 
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question whether nuclear forces are, partly, many particle 
forces. 

(2) They may help to decide between certain ambiguities 
( such as the exact space and spin dependence of nuclear 
forces ), which the data at present available do not settle. 

From the very fact that more than two intimately inter- 
acting particles are involved it ts clear that collisions of 
composite nuclei present a very complicated theoretical problem. 
While it is encouraging that the corresponding bound state 
problems have been solved with a reasonable amount of effort to 
within a few per cent, it may not be out of place to mention 
the additional complications of the scattering problems. The 
asymptotic form of tne wave function, instead of having to 
vanish at inflnity, must correctly describe the colliding 
particles. Above certain energy threshholds it must exhibit 
the possibility that part of the kinetic energy of collision 
may be used to raise one or both of the colliding particles to 
an excited state; or that an exchange of nucleons may take 
place, so that the resulting nuclei differ from the colliding 
ones; (@ef- 1) es) We? 4! )3 or, finally, that a triple 
or multiple disintegration may take place (e.g. ae ae 
( intermediate stages ) —> AHet+n' . )- Interactions 
with the radiation= and electron-neutrino fields will not be 
considered in this paper. We shall also exclude from our dis- 
cussion proceades in which multiple disintegrations are 
energetically possible since in that case the asymptotic form 


of the wave function is very complicated and it was considered 





advisable to concentrate first on a solution of the simpler 


problems. 


- Experimental Datae 


In the last twenty-five years a great many experimental 
investigations of collisions of light nuclei have been carried 
ae Many of those performed before 1937 are discussed in aie 
Bethe's and Livingston's review articles on "Nuclear Dymanics" : ° 
A few other typical references are given Pick. Scattering 
experiments with neutrons as of 1947 have been compiled by 
He He Goldsmith et + Quite recently, collisions at very 
high energies have been investigated by means of the Berkeley 
Sveloiwen 

It is our gerie ved aim to make possible a theoretical 
interpretation of these experiments. However, this paper is 
devoted mainly to the development cf theoretical methods and 
only very few numerical calculations have so far been carried 


out. Hence we shall not be concerned here with a detailed 


analysis of the experimental data. 
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10. He Ae Bethe, R6Ve MOS. Phys. on 69 (1937). 
ll. Me S. Livingston and He. A. Bethe, Reve Mod. Phys. 9, 245 (1937). 
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Previous Theoretical Work. 


Theoretical work on more complicated scattering problems 
has been mostly of a qualitative and phenomenological nature. 
It has been largely based on the resonance formula of Breit 
and bainn | oti exhibits the scattering process as the 
temporary formation of a compound nucleus with subsequent 
emission of the incident or some other particle. In terms of 
this theory it has been possible to understand the general 
features of nuclear scattering, but no detailed information 
about nuclear’ rorees can be derived from it. 

ntien ore Heisenberg Uneve ave yzen the general 
framework of collision processes without going into any detail 
of the interaction. Their theory serves two purposes: It 
elucidates the general roatneee of scattering phenomena, 
regardless of the nature of the interaction, and it represents 
the first stage of any detailed calculation of composite 
collisions, based on some specific nuclear theory. 

Furthermore, as an approximation to the correct wave 
function, Wieeler "ia used functions of so-called group 
structure, which have been employed in all succeeding 
calculations. To determine the best possible functions of 
this sort, he made use of a variational principle. The 
contents of sections on and I.7. are an extension of this 


methode 





17. G. Breit and Ee. Wigner, Phys. Rev. 49, 519 (1936). 
18. Je Ae Wheeler, Phys. Rev. 52, 1107 (1937). 
19. W. Heisenberg, Zs. f. Phys. 120, 513 (1942). 





More directly related to the present paper are the works 
20 aN LU 23 
of Schiff , Ochiai , Fltgge , Motz and Schwinger , 
2427S” 2b 
Buckingham and Massey » and H8cker who have calculated 
scattering cross sections in n-d and p-d collisions from 
approximate solutions of the Schroedinger equation. 


These calculations are all based on a wave function 


with group structure, that is a function of the form 
(O.1) NO y : V5 Fi: 
where Y: » Ys, describe two collidors 1 and j when 
outside their range of interaction, and Fi describes 
their relative motion. The sum is taken over all pairs of 
nuclei whose formation is compatible with the given energy 
and over all permutations of identical particles, with the 
signs adjusted in accordance with the exclusion principle. 
The functions fs are then chosen so that ¥ satisfies the 
Schroedinger equation as well as possible. A function of this 
type has the correct symmetry properties and asymptotic . 
behaviour; but in view of the strength of nuclear forces one 
cannot expect it to approximate too well the correct _— function 
in that region of the configuration space where the collidors 
Interact, and it is just this region which determines the 
scattering cross section. 
20. Le Ie Schiff, Phys. Rev. 52, 149 (19357). 
21. Ke Ochiai, Phys. Rev. 52, 1221 (1937). 
ele S. Flugge, Zs. f. Phys. 108, 545 (1938). 
23- Le Motz and J. Schwinger, Phys. Rev. 58, 26 (1940). 
24-6 Re Ae Buckingham and He. Se W. Massey, Pe. Re So, 123, 

179 (1941). 


25¢ Re As Buckingham and He S. We. Massey, Phys. Rev. 73,260 (1948). 
C66 K. H. HOcker, Phys e ZSe5 43, 226 (1942). 
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Of the calculation mentioned above, those of Buckingham 
and Massey ee the most extensive. Their method is in 
essence the following: 

The centre of gravity motion of the two neutrons and 
proton, participating in an n-d collision, is separated out 
and there remains the following equation for the wave function 


in the relative coordinates: 


(0.2) ATs Vin. + “234 Vail ¥ = {E,, T E a} y 


Here T is the kinetic energy operator for the relative motion 
of the eines particles; Vig i3 the interaction potential 
aterogsa particles i and j, supposed independent of the position 
of tne third particle; Ee is the energy of relative motion of 
the colliding neutron with respect to the deuteron in the system 
in which the centre of gravity of all three particles is at rest; 
and Es is the binding energy of the deuteron. 

Let us introduce relative coordinates to which we shall 
have frequent reference in the following pages. If chy oe 
and R, are the absolute position vectors of the two neutrons 


and the proton, then C, » 8a, 93 are defined by 


Sis R.-R, 
~- 9. = RRs 


3 = R,- Re 


and similarly Y, ff, Y, are defined by 





“2 
ma 
— 
— 
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RitR; 
"= R, 2. 
(4) v2 R, — Ryehs | 
ze 
R, +R 
= — a 





Fige 1 


In terms of these coordinates T has the form 


Ue 
K wv —2L 
T = val a Z f, 
or the two other forms obtained by replacing 1 by 2 or Se 


Mis the mass of a nucleone 


as 


The wave function ys is a function of the six relative 
coordinates represented by F and r ( or any two linear 
combinations of these ) and of the spins of the three particlese 


Equation (0.2) is to be solved subject to the condition 


— ss 


that when either Y, or n become large, eo describes a 


deuteron and a free neutron colliding with ite Thus 
ot 


(0.6) y PA: o(¢,) ie + —— )x, ae r ee 


™ lm) (et" + hg Xe & ¥, > co. 
| | 
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where <9 (¢) is the ground-state wave function of the deuteron, 
he |k\ = (s%" Ey funy , y= VR, = | 
scattering amplitude and i239 ‘Ya. are appropriate spin 
functions which differ according to whether we consider quartet 
or doublet scattering. 

Although we are not familiar with a rigourous proof that 
a Six-dimensional partial differential equation of elliptic 
type has a solution which is finite everywhere and satisfies 
such asymptotic boundary conditions, one intuitively expects 
this to be the case. However a reasonably rigourous solution 
is undoubtedly very complicated. Buckingham and Massey, to 
simplify their problem therefore assume that the exact solution 


can be everywhere well approximated by a function of the form 


7 © = ag) #(n) q+ ¢(e) €(n) He- 


This is equivalent to allowing only for the distortion of the 
free neutron function ( since F is as yet undetermined ), but 
disregarding the distortion of the deuteron by the neutron. 

To find which function F makes VY the best possible 
solution, they substitute (0.7) into (0.2) and, using the wave 
equation for the deuteron, they find that F should satisfy, in 


the case of doublet scattering, 


(awe + Ey - o- Div, e Mn) \ F Ln) ¢ls) 


(0.8) 


: c= SANs, al9s) FLY) + Va3 ¢(.) F(a) 





Here w, h, m and b represent the strengths of the Wigner, 
Heisenberg, Majorana and Bartlett forces. 

For wae fixed value of oy » (0.8) could in principle be 
solved for F (rv) giving Flr). But as there is no reason why F 
should be independent of g it 1s impossible to satisfy (0.8) 
by a single function F(r). This is just an expression of the 
fact that a function of the type (0.7) cannot solve the wave 
equation exactly. 

To obtain, as it’ were, an average solution for F, 
Buckingham and Massey multiply (0.8) by (9,) and integrate 
over the 9 ~-spacee This leads to an integro-differential 


equation for F, of the form 


(0.9) (WH) Flr) = (2v-h) UC) F(r) | 
+ (Zu ) { Q (v, r'} F(v')elr'), 


where :U and Q involve the known interaction potentials and 
the deuteron function. This equation is finally solved by 
numerical methods. 

Previous to Buckingham and Massey, Motz and Stieingeh | 
have calculated the scattering of neutrons by deuterons in 
the zero energy limit. The essential features of their 
calculation are similar to those of the British workers. 

Thus, they too consider the deuteron as rigid and are led 
to an integro-differential equation for the function F. 
To solve it they replace it by an integral equation according 


to the schemes: 





10 
(44%) Flr) = ( K(y, mJ er" )(Ae') —_ = 


(0.10) , Pill ae tf 4 \ ty ( 
F(r) = cael \ 3 Kr Ela el (ar), 


After carrying out angular integrations analytically the 
remaining integral equation is solved approximately by 
replacing it by a set of simultaneous linear equations. 

painaliy, Beker” has tried to write Flr) as @ power 
series in*wand determine the coefficients to best fit the 
Schrosdinger equation. However his results depend, by a 
factor of about 1.5-2.0 on the point where the power-series 
is taken to pass over into the asymptotic form of F. 

The other calculations do not contain any essentially 
different methods. 

The different authors using various kinds of ( central ) 
interactions obtained widely different results. For the 
total cross section of ned scattering at O-energy they are in 


10 z{ 
units of barns as follows: Schiff * 6.11 3. Ccheal Gao. 22: 


Motz and Somme sraeees 4.57, 6.91 $; Buckingham and acco 
5eO, 3.9 ( extrapolated ) ; donne os 1.88, 2.95 . The 
experimental value is near Se 

Other two- and three-particle systems do not exhibit such 
a strong dependence on the type siplaeenccie one at low energies. 
If taken at their face value, then, these results aula 
provide a useful method for distinguishing between various 


kinds of interaction, all of which are compatible with other 





27. Le Je Rainwater et al., Physe Rev. 73, 733 (1948). 





al 


low energy phenomenae 
However, sucn a procedure would be without firm foundation 
in view of the unreliability of the group structure function, 
on which all calculations are based. This emphasises the 
need for a systematic method of solution, which is not 


restricted to any particular kind of function. 


Summary of the Present Paper. 


_It is not the primary purpose of this paper to add 
another approximate calculation to those already published, 
but rather to investigate the possibility of a systematic 
procedure for solving the Schroedinger equations corresponding 
to composite scattering processes. It was natural to look 
for a variational method, as this has been the chief tool 


for solving the problems of the bound S- and 4-particle 


U-b, 28 2.4 
systems and has recently been applied by Schwinger 


to 2-particle scattering problems. 
There are two types of practical variation procedures: 
One is the Rayleigh-Ritz dietheaen oe undetermined coefficients; 
the other the variation iteration iat eom. one Green's functions. 
We have tried both methods in the many particle scattering 


problem, and we strongly incline to the opinion that the 





28. Ne Svartholm, Thesis, Lund, Hakan Ohlssons, Boktryckeri (1945). 
29. Unpublished lectures. 
30. We Ritz, Je reine angew. Math. 135 (1909). 
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Rayleigh-Ritz method 1s the more useful. 

Using this method we have made exploratory calculations 
on 2- and 3-particle scattering problems, some with the 
inclusion of tensor forces. On the basis of these we believe 
that this procedure can indeed be used for a systematic 


investigation of collisions of light nuclei. 
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ft. THE RAYLEIGH-RITZ METHOD IN COLLISION PROBLEMS. 


Iele General Remarks. 


In his original paper on wave mechanics, Schroedinger 
derived the wave equation from a variational principle 
analogous to Hamilton's Principle in classical mechanics. 

The variational formulation of the equations of motion of 
wave-fields has been very fruitful not only in theoretical 
developments of field ewe om but also in approximate 
solutions of complicated quantum mechanical prone roa 

The method most extensively used in practical problems 
has been the Rayleigh-Ritz eee: linear trial functions. 
However its application has been limited to systems in bound 
states, whose wave functions vanish rapidly at infinity. 

Reeently, Heteenteeg™ has pointed out the correspondence 
between the phase shift energy relation of a collision problem 
and the energy spectrum obtained if the participating particles 
are enclosed in a large box. 

Consider, as the simplest example, a particle of mass M 
under the action of a short range scattering centre. The wave 


equation for the radial part of the spherically symmetric 


solution is 





Sle Ee Schroedinger, Anne Ge Phys. 79, 361 (1926). 

32. See, for example, G. Wentzel, “Einfllhrung in die 
Quantentheorie 4. Wellenfelder", Wien, Franz Deuticke 
(1943). 

S5e Ee Hylleraas, Z3e f. Phys. 54, 547 (1929). 
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(1.85 \ i £ + w (>) ula) = Eula) 


where U (n) is the scattering potential and E the total energy. 
Asymptotically the solution becomes 
ea 

ec} te ——= cin ame p+ 4 (€) ) 

h* J 
where the phase shift n(e) » corresponding to the energy E 
of the incident particle must be found by a detailed 
Solution of (1.1). 

Now let us ask for the energy levels of the same 
particle in the same potential field when enclosed in a 
spherical box of large radius, a, on which the wave function 
must vanish. Since the differential equation is unchanged, 


all those solutions of the scattering problem will satisfy 


the bound state problem for which 
1.6 Jane. ( <) = 
( ) | S Ww ( Fe a + 1 i O. 


It is clear that if we have solved the scattering 
problem and determined the func tion 4 (€) we can find the 
energy levels, Ee. of the bound particle oe, solving (1.3) 
for E3; and conversely if the spectrum, Es of the bound 
particle has been determined the phase shifts y (E;) can 
be found from a GD 

| These considerations lead us to expect that the 


Rayleigh-Ritz method, which has been so valuable for 
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calculations of energy levels of bound systems, will apply, 
| with certain modifications to the calculation of phase shifts 


and scattering amplitudes. 


Ie2. Three Related Treatments of One Dimensional Problems. 


The Collision Determinant. Numerical Illustrations. 


In this section we shall illustrate our methods by 
treating the scattering of a particle by a fixed, short 
eves scattering centre. We shall show that this problem 
can be attacked in three closely related ways. 

Let us first simplify the Schroedinger equation (1.1) 


by introducing a convenient unit of length, A, , and writing 





“ 
(2.1) x= 7 ) ee aot De : V (x) = a aM mt a 


It then becomes 


(are) |S e vx) | te 


dar 


It will be our problem to solve (1.5) subject to the conditions 
that 


(25) ulo)zo, u(x) —2> Sin kx 4a. y Cor KX 


B-y “C a <—O e 
In particular we are interested in the phase shift which 


determines the scattering cross section. 
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Method 1: The Energy as Proper Valuee 

We know that (2.2) and (2.3) have a solution for every 
positive K .« However, instead of looking for the phase shift 
corresponding to a given K , we can reverse the problem and 
ask for the K corresponding to a given logarithmic 
derivative, L, of u at some large distance ae This means we 
look for tne proper value is belonging to the differential 


equation (2.2) and the boundary conditions 
(24) ulo)2o , wla)=b u(@). 


34 
This is a proper value problem of a standard type and it 
= 
is easily verified that \K as calculated from the equation 


(225) ( \ ui + Kae Ver | ey L u>(a) “ye 
0 


is stationary with respect to all variations of u which 
satisfy the condition 5u(0) 2U. 


The Rayleigh-Ritz method is directly applicable to (25). 


We simply try for ua linear combination of trial functions 


(206) w= u(x) +4 Spy VID ics oe Cu Uy, (x) 


( with all LAS (0) =0 ) and substitute this expression into 


(2.5)e This leads to an equation of the form. 
2 ~? Vv 
(27) Qed (Ai, +82, 8 ) e: Cc) =0 
az 
y= 


where h:. Pi, are integrals involving uz ond ONE 
Pei Se TR ne oe 


34¢ Re Courant and De. Hilbert, Methoden de Mathe Physik, 
vole II, 508, Borlin, Springor, 1957. 
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Lv 

The stationary property of WK with respect to variations 

of the c+ leads to the following system of equations for the 
determination of these coefficients: | 


* an : 
7% 22 (Ay + Bi; k Jeseo Aelydy.-en. 
a" oC, } 


As in the case of the bound state problem this system is 
compatible only if wo has certain proper values Ko which 
can be found by the solution of the corresponding secular 
equation. They are approximations to the correct proper values 
of (2.2) under the boundary conditions (2.4). Their errors 
are of the order of the square of the errors of the trial 
functions. | 

Now the phase shifts of the exact solutions corresponding 


to the given L satisfy the equation 
(2.9) Kk cot{ Ka+n(w)) =b 


If we substitute our approximations Ke in this equation we 
obtain, to the same (' namely the second ) order of accuracy, 
values for a] Ky) corresponding to the energies Ke 7 
Then, either by varying the pre-assigned value of L or the 
"cut-off" distance a, we can obtain a for all other energies. 
In this form of the variational method, we have followed 
very closely the analogy between bound and unbound systems 


pointed out by Heisenberg. However the quantization of the 
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energy which entails the necessity for solving a secular 


equation of the usual kind seems artificial and can indeed 


be avoidede 


Method 250 Ths Logarithmic Derivative as Extremume 

| Since for a given L, > as calculated from (25) is 
stationary, so also must L be stationary when calculated with 
a given a Using the trial functions (2.6) we are led 
again to (2.7), but now we take K as fixed and L as 
adjustable so as to make the ensuing system of equations 
compatible, Therefore we exhibit L explicitely and write 


Q in the form - 


Qe = Cz, UK Jet L , C,4,(a) +c, u, (9) 


(2.10) ‘y= 
v 
oe CL Ua (a) { =e) 
It is convenient to transform the coefficients c- to a new 
set da a as follows: 


a Cm, (a) 48 C ee, Eee eee CY uy ( a) 
(2.11) 


diaz ec: | VY 
In terms of the d;, Q has the form 


Az\ 


(2012) Q = Z us 3 d A\ +L d =— 70 


y= 
The stationary property of L with respect to variations of the 


ds now leads. to the equations 
A 





i? 


Zz Dy (er) dy +b dy =? 


y=! 
(eels) 


a Dry (x) ok} =O 
joi . 


= 


) AS Ly 3, <2 ee WY. 


To make the system of equations (2.13) compatible L must 
satisfy the determinantal equation 


~~» ey as 


hee 4) ie 


e —_ om @& 
. : 





Hence L is a quotient of two determinants, each of which is 
a function of "ke 


Let us note that A is symmetrical and of degree l 


in L, in contra=-distinction to the secular determinants in prob=- 
lems of the bound state, whose degree in the unknown energy 


is as high as the number of trial functions. We shall see 


later ( sections I.5e ~ I.7. ) that the structure or A or 


of the corresponding determinant occurring in method 3 reflects 


the type of the collision and for this reason call it a 


3° 
collision determinant e Its degree in L is of course 


oe ee 
oSe It is related to Wheeler's Fredholm determinant, ref. 
fe ie, eqe (45). 







—— = 







= 





>_> S_ =-> = > 
a @ @ Gees «: 
<< & - 
om 
S_ i waa 
2 a “ae | 


aS ——= 
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invariant under the transformation (2.11) which was introduced 
only for reasons of greater clarity. Indeed one can easily 
verify that the determinantal equation for L, as deduced 
directly from (2.10), 


i835) Cyn (a) wy (o) | = 


iEmeaiso of the first degree in 0, since cre arr 
\\ wi (a) Uy (a) \\ of the coefficients of L has rank le 
There remains as an unsatisfactory feature of this 
treatment the dependence of the phase shift, deduced from 
ee and. L, on the cuteoff distance a, even when a is 
well beyond the range of interaction. This question will 
be discussed in more detail further on in this section. 
feanwhile we shall indicate another method which does not 


have this drawbacke 


Method 3. The Tangent of the Phase Shift as Extremume 
Let us return to equation (2.2) and ( writing € 


mirnstead of K~ ) note that, in the bounds gate,ae 


calculated from 


(eac) TL =z \ « “fT ae V(x) } 4 da = 


is stationary relative to all variations of u with §u(o) => 


Su (co) - 0 -« For, 


(2.17) 





oe 


and both the integrand and integrated terms vanish. 
In a collision, however, the wave function has the form 

(2.3) at infinity. If we restrict our trial functions to 

those having the correct asymptotic form, although 


~ 


possibly wrong phase shift, then 


(2018) Su —> Stan . w> KX an x —> co. 


Hence 


A ol 
(2019) LA ju Se Su Fu | 


7H= 
so that 


(2.20) § e=: KK Ann 4) =U. 
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Since for the correct wave function I =0 » it follows 


that the equation 


i2e21) T+ k Aon 44 = K AQwy | 


where x is the trial phase shift, is a stationary 
expressicn for tan 1 ° 

To apply the Rayleigh=- Ritz method it is essential to 
write this variational principle in a homogeneous forme 
Let us therefore consider wave functions of the asymptotic 


form 





S60 Le Hulthen ( Extrait, Dixieme Congres des Mathematiciens 
Scandinaves, Copenhague, 1946 ) has used an equivalent 
equation to find 4 by further restricting his wave 
functions to satisfy the condition I=0. 





Le 


(9.02) wu ——» A sim Rx + Bur ¥x 


instead of (2.3). We then find in the same way as above 
ie 
(2.23) T+ ABrs BD K tay 


as a homogeneous variational principle for tan " - If we use 
a trial function of the form 


2.24 = 
( ) UL AswwHkRwt 4 C, FEE. 


where 


A 


then since ot multiplies only ee » Substitution of (2.24) 
into (2.23) will give a quadratic form of: the same structure 
as (2.12) which in turn will again lead to a linear collision 
determinant for tan ° 

The drawback of this formulation is that the trial 
functions must be of the form (2.24), (2.25). For practical 
purposes it may therefore be preferable to re~adopt a finite 
cut-off as follows: Let a be a distance beyond which the 


interaction is negligible and let us use as trial function 


Ae = ae = Cyp4, + CC, bh, 4+ fee CO Un ) ae 
(2s 26) | . 7 
w= rol = Bh SAnKkx a Bes kx 
where a 
(4 V(t) 
A= { «'(a) age + % A [a). coxa | 
(2.27) 


| i , 
a yw (9) wt KQ ~ uw’ (a) Sim vay 





ZO 


A and B are defined so that u is smooth across x =a. Clearly, 
Since u has the log, Ge xX 7 a, the intesrand. or l Bmore 
for x pat. ( see (2.16) ) so that 

Saye at oy 
(2.28) re Je + e- viv) ts ot . 

o 
When (2.27) and (2.28) are substituted in (2.23), we have a 
variational principle involving only uw? (x), whose asymototic 


form no longer concerns us. The price paild, of course, is the 


reappearance of the cut-off distance a. 


The Well-Depth as Proper Value.* 


Before going on to a discussion of these methods let us 
note that if we write V(x) =< -h\ P(x) » where f£(x) determines 
the shape and (the depth of the potential well, we can 


regard 1‘ as the proper value in the Schroedinger equation 
ue 
(2.29) co re + “he f(x)} u (+) = 08 
dx | 


A) as calculated from the equation, 
~P 


(2.30) ( u (x) 4 a " ame aN 4 (x) \ u(x) 20 


is stationary if 


(2.351) $u(0) =O ) Sul ~ } oO 





% This tdea was used in connection with bound state problems 
by Svartholm ( ref. 28 )-. 
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Both K and 4 are now regarded as known while the strength of 
the interaction is the unknown. This is in contrast to 
methods 1-3 in which the interaction was regarded as given 
and and 4 or K and L were determined as functions of one 
another. 

While this approach has no advantages in connection 
with the Rayleigh-Ritz method, it provides a useful view- 
yvoint in the variation iteration method to be discussed 


later ( section II.8.). 
Connection vetween Methods 1-3: 


isthods 1 and 2 differ only itn respect to which of 
the quantities L and Ko is presumed known and which is 
ericnown. With the same trial function and the same a they 
must lead to the same relation between tank andl Ke. 

On the other hand method 3 differs slightly from 1 and 
oe Let us trace the connection between the cut-off independent 
form of method 5 and method 23 

In the latter, lL is determined from (2.5), which after 


an integration by parts gives 
ghas - 2 
ss) - uz, | . tT +b wis) = 0 
x= 


where I is defined by(2.28). If u has the asymptotic form 


(2.22)we can evaluate the two boundary terms, which leads to 
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-—— sim (Ka +4) m (Ka+ 4) Fak 


Co” 4A 
+ —— 2 cor | Kat Siu Kat") =D 
is "1 "+) “a 
where "x= tan ( 8/4) e Multiplying (2, oo) by ‘ *(Ka+y ) 
° pis a 


we obtain 


(2.34) 


px. | cA (Kat 4) - UA [Ka +ye)f + 
: Py quasde Lm =a 


Sib (Kats) 
At the particular cut-off distances where 


(2.35) Ka = ~htl ay 
2 ) 


(2.34) becomes 
(2.36) a. + Aw Tow "se = A’ 2 tan Wks 


This agrees precisely with tke result (2.23) of method 3. 

We now ask how h as determined by method 2 varies with 
a for a given trial function u which has the correct asymptotic 
form but incorrect phase shift. In (2.34), let us use the 


abbreviations 
4, = Kot Le dy = qa 


(2.37) = 
J = = we 


By some elementary manipulations we can then reduce it to 
(2.38) WX $y = oA 4 = ay 


For the correct wave function, J =~ 0, 30 that by (2.38) 





Kins 





r 
ais € 299 = hg 905 «6') HAYES 


2°O 


¢ 


roan 





26 


Sy zo » except where cot 2 = + oo, that is where 
KAM, = nw « At those points $y is undefined. This 
is understandable from (2.5) since at those points u = O so 
that Lis me determined. 

The variation of Jy with z is shown for different 
values of J in Graph l. For small values of J, cot 2 is 
negligible as compared to 1/J, except at the danger points 


KQ4ne=9 = nW » 80 thet in the limit as J —90, 54 


becomes independent of a. 
Bilinear Formulation. 


The Schroedinger variational principle for bound systems 


is usually stated in either of the two forms: 


Q \" f- (SE) - ve) ats ew) ola 20 


cil 


(2.39) | wat 
SQ=0 
or - 9b Ae 
= hats le dx = 0 
® @2), 44s, lx) +e) « 
(2.40) 


SQzO 
where all trial functions satisfy 
(2.41) du(e)=o, du( xe) - 0- 


Let us recall that a modification of (2.39) for the 
case of a collision led to methods 1 and 2, while (2.40) 


was the starting point for mothod 35- 








ae 


We can replace (2.39) and (2.40) by the following 
equivalent statements, involving bilinear instead of 


quadratic forms: 
oD 


Q, = am bey dur ol @ u, +t Ecu,u | dx 20 
(2.42) r 9 in) a) Ux) . ae 

dQ, <0 
or Db 

d* | 

| u | < ~via ef dx = D0 
(2.43) ° ea 

§Qiz =v 


where all trial functions U, and U, satisfy 


\ 


(2.44) duzlo)20,  §ule) = o are 


* 


) 
3 in (2.42) and (2.43) indicates independent variations of 
Uy and U,- The truth of (2.42) and (2.43) is directly 
verified from the Schroedinger equation and the boundary 
conditions (2.44). a | 
Clearly the corresponding variational principles for 
Peliision problems can be similarly modified. Without 


going into details we shall merely state the results. 


Methods 1 and 23 


se SCAM) 2a V0) 


+L u,(a) ur (4) = 0 


Aria se 
if the trial functions satisfy 


(2.45) 
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(2.46) §u:(2) 29 4 Ae Wyte 


CO Compames to (2.5).) 


Method 33 


ee = AF Sys aA ~ A, Aron y 


whut ee J \ = — 4K ~ V(x) ft 


(2.47) 


oly ® 


(2.48) 8 Kee 


if the trial funetiongs satisty 


du: (0) = | 


(2.49) A= ee 


UU; —— > Az sinnx + Bice ka 
( Compare to (2.23), (2.24).) It is not at once apparent 

that tan 4 as calculated from (2.47) is symmetrical in UU, 
and Ws That this is so, however, follows from the fact that 


(u J2 xt v(s} uy de 


° 


ey 


Fu + te Vl (x) \ ay te (ony Zz), 


ra 


(2.50) © 


c( 


Lin + A, 6 a Ps, 8, Kk. 

In one dimensional problems there is no advantage in 
using the Rayleigh-Ritz method in conjunction with the 
bilinear variational principles instead of the quadratio 


Ones. But we shall find the bilinear version useful in three 
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dimensional problems ( sections I.4. - I.5. ), where the 
Schroedinger equation has more than one solution, and in 
connection with the variation iteration method ( sections 


ie Oe a Li sk. es 


Numerical Tllustrations. 

Using the methods 2 and 3 we have calculated 23 
scattering of neutrons by protons at zero energy. We have 
assumed the interaction 


V(x) = -b.or4as Si 
fe, ) 


= ao xX>] 
( x measured in units Of Ay=2-30 ian em ), which gives 


the correct binding energy for the deuteron. The wave 
equation, in this case, is simply 

yr 

Au 
(2.52) (= 4 Vla)u =O, 

ee 

and can, of course, be solved exactly. For x > 1 the 
solution has the form u =A (x + X ) where X is the desired 
limiting ratio of tan y /*&, 88 K approaches zero. 


As trial functions we have used 
h 


(2.53) uy =z 2 Cyx* 
! 


with n ol, 2and 3. Since (2.53) does not have the correct 

asymptotic form, it was necessary to use method 3 with a cut- 
off, a. The interaction vanishes at x=], which is the most 
appropriate choice for ae However to study the dependence of 


X on the cut-off distance, we have made calculations also for 
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larger values of a. Clearly a function of the form (2.53) 

can approximate the correct wave function better in a small 
than in a large interval, so that as a increases we expect 

our approximation for X to get worse. 


The results are listed in Table I: 


Table I. 


Approximations to X & Aires ton 4 [\« ‘ 





| > he 3 
| estest faxtest tes 
Yethod 2] Method 3 


- 2.084 














+ 1.3558 






te Le 126 






+ ae SOO = 4.039 







+ 1.358 - 51.606 





+ 1.3538 + 7.654 


These results exhibit a number of significant features: 


1. The methods 2 and 3 lead to errors of the same 
order of magnitude. ( In the present example method 2 is 


somewhat better. ) 


Ze As expected, the results grow progressively worse 


with increasing a. 


Se The dependence on a becomes weaker the larger the 
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number of parameters ( except for the 1.338 column; see 5. 
below ). 

4. There appear singularities, owing to the wrong 
asymptotic behaviour of the trial functions. As the number 
of parameters increases the singularities move outwards, to 
Poa aaity. 

5. The trial function CX; in method 3, leads to a 
result independent of a but of the wrong sign. This can be 
understood as follows: Ce x is the correct solution corresponding 
to no interaction. It has ( accidentally ) the eorneee 


asymptotic form ( though wrong X ) and therefore, by the theory 


of method 3, must lead to results independent of a. In fact 


we have a 
LY 
aaa | - V(x) Uy olic 
(2.54) a 
\ = 
B = ma C1) 3 = OD 


all.of whicn are independent of ae Using these in (2.23) we 


“obtain for ee determination of X = tan "IR 4, 
(2.55) ~ ( Vr (x) tx = ¢,” tan y /K 


which is precisely the Born approximation for tan . For 
small V the latter gives good results. But as V increases 
( tan ” Bon, has a fixed sign while the true tan " 
eventually changes sign. This has occurred in our case, 
hence the wrong sign in the 1.558 - column. 

In summary we may say: 


(a) Method 1 is the least natural to collision problems. 
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(b) Method 2 appears to give slightly better results 
Bunan method 33 but it suftersmerune aco flaw of necessarily 
involving a cut-off distance. 

(c) Nethod 3 is closely related to method 2. It does 
not necessarily involve a cut-off and yields the Born 
approximation for the phase shift as a special case. 

| (d) The Rayleigh-Ritz principle applied to either 
method 2 or 3 leads to a collision determinant instead of the 


usual secular determinant. 


I.d. Variation of the Phase Shift with Energy. A Stationary 


Expression for the Width of a Resonance Level. 


We have seen that methods 2 and 2 lead to stationary 
expressions for the phase shifts corresponding to arbitrary 
energies. It follows that it must also be possible to 
obtain a stationary expression Tor ay [ A K which determines 
the variation of the scattering cross section with energy. 

It is indifferent which of the two methods we use, but 


method 2 seems more convenient. We introduce the definition 


Q. 
en Me(a) 2 J \- Fe) + ee vie sea), 


where Ly is the logarithmic derivative of the correct 
ee 
function YU, belonging to the energy Kk . By the Schroedinger 


equation (262), My (uy) 20> and by the stationary property, 
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(2.5), we have rigourously 


(3.2) My lu) = My (Sux) = 0( Sux) , 


mor any function U=Uy tduy - Differentiating (3.2) with 
respect to we obtain, to the first order, in Ou,, 


EB sy * Ma(uy O 
A 14 


Hence by (3.1), we find the equation 





Q 


ie 
(3.4) 2x | n ge 4 Fe u (a) =0 





0 


wolch is a stationary expression for Alby JA 
z 
From this equation we can determine dy /a(ui)as follows. 


We note that 
(3.5) Ly = Ker (Katyn) 


Hence 


v cli 
ey, ayy = et AY lm) ee (Ka+ya(ua +x =I) 


so that 
dun es a “ K 
A>) = AK 
(3.7) ALK _ Katy x) , taleary a) _ 
7 be Aes u Ke 
q 
| ar lx 
= ts sit (da thn) 4 sink tlea+nell 
uw u(a) ” Sf ae 


a 


o 
— ——— ° 


2 







i Te 


oe ee, - 









-— 
oom — 


o4 


We may observe that this stationary expressrong er 
involves only the wave function and the correct phase shift 
Wk ; ale sae reference to the interaction has disappeared. 

Let us check that if dy, [dld)is computed from (3.7), 
using the correct wave function U,, the apparent dependence 
of dy, {alx’) on a disappears, as it must, as long as a exceeds 

0 


the range of interaction, Q,. 


For x > Q,y , the correct WY, has the form 


(5.8) Uy(x) = C Sin(Ke +4) 





Hence 
Q 
(3.9) Se eS ak + 
dle) wer? eee Qk 


Differentiating with respect to a we find 


cea d(K)) ~ Fou Sin (a tn) + = eat ¥, 


(3.10) 


cl 


Lun afen mal ate2teney) 4 
ae 
which proves that dy/d(w) 1s independent of a. 
AS an application of (3.7) we shall derive a stationary 
expression for the level width at resonance. 
In the nelghbourhood of a resonance energy, defined by 
the condition that 1K = As) T » vari "ie has a singularity 


of the form 





5S 


Vem 
(3.11) Ton 4 4 7 aan 


v x +) 
Ko- Ky 


L 
where |<, is the resonance energy. {' 4s called the level 
width because the cross section in the vicinity of K, is 


given by 


w (GT) 


5 12) 3 oO = ¥ . 
meen) eS) 





4 


Now near Kany 
of 
(3.13) aks = [na s\n ~b = (K = Ky) , 
| A (x,") 
Therefore 
| Ste 
(3.14) Tan 4 x > = Be: | . 


L 
Ko KY 


Hence, comparing (3.11) and (3614) and using (3-9) with y= laeS ) 


we find the following stationary expression for 4 " $ 


t 











are 
“is 
fis.15) AC KS) 
; UK 
= 4 Ka 2 : Lv | ze 
ua) Ka. | ol < AAK ~ Sin(A aX) 
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1.4. A Variational Principle for the Scattering Amplitude, 


When dealing with a wave function of given angular 
momentum, a correspondence to a bound state problem can be 
established because we know, a priori, that asymptotically 
the solution has a radial logaritnmic derivative independent 
of the angle. But if we want to treat the whole incident 
plane wave at once no such information is available. We 
have already seen that only method 3 need not involve an 
enclosure and that the Born approximation forsindi vidual 
phases can be deduced from it. Both facts suggest that it 
“would be the most promising attack on the scattering amplitudes 
of three dimensional problems. 

The wave equation for the relative motion of two colliding 


particles of equal mass M is 
kg 1 = 
aa f-S 9+ U(R)} LIR) = = YL (R), 


where R is their relative position vector and U (R) is the 


interaction potential. Setting 


R Z fre » 
re Re BE yy Mle. 
42 





(4.2) 
Wr) = (ao¥) 


where Ny, is our unit of length, it becomes 





oT 


(4.3) Lv" ee v(e){ (+) =D 


| We shall now follow a procedure corresponding to the bilinear 
form of method 3 ( see (2.47) ~ (2.49) ). We define teu, » and 
Kp as the two solutions of (4.3) which correspond to plane 


waves incident along -K, and \<, and have the asymptotic forms 





ake. he oe 

fei, ae eg {omyn) SO 
(4.4) : 2 
Re en" 

Soe >e vr +t {[¥, gk) = 


Here (Ky \ = [%| =[¢) and { (2, %) is the scattering amplitude 
in the direction IK of a wave incident along It: . Corresponding 


A 


to (2.47) we now define 


-—) 


(4.5) ai 7 : ; | y. i, & si i v) 1 kk, (ar) 


d 


In virtue of (4.3), we find after an integration by parts, 


ao) St Ky, K, = \ Wise, Ihe, a SK, = ¥-«,| as 
S 


where S is a large sphere, and 2 denotes differentiation 
along the outward normal. 

Now let us admit only trial functions of the asymptotic 
form (4.4), but with possibly incorrect f's, so that 
asymptotically 

aKA 
(4.7) Seu, il 54° (%, «) ae 
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| (vt) 
where 3 f (K,, i) Ls the error in the scattering 
amplitude of the trial function Yu, « Clearly only the plane 


wave part of Vu, contributes to the surface integral (4.6) 


so that 
Avr 
| oth cn | 
$1.4, > er 
(4.8) S | : - . L : 
| = ave a on TSE, K) dS 
A Orn 


3] () 
Dirac has shown that the factor multiplying jf (K,&) has 
a § - function like character for large values of N, s0 


that the integration gives simply 
_ ¢ g@) 
(4.9) § Liu,, 4, => Ui 3 { CAE 


memice Yor the correct Yr, » age ete O » it follows 


that f(k,, 4) » as calculated from the equation 
(rv) ‘ ~ 
(4,10) am Kw, UT 4 Gt OG Ki} = KW 4 (k,, Ky) 


is stationary with respect to independent variations of YW -K, 


and Yeu, ? subject to the conditions (4.7). 


Symmetry Properties. 
Equation (4.10) involves two trial function, TK, and 


Ya, - These differ by a rotation in the coordinate space 





S7e Pe Ae Me Dirac, The Principles of Quantum Mechanics, 
( Clarendon Press, Oxford 1947 ), third ed., p. 191. 
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and, in general, also in their functional form ( to the 

first order ). Let us indicate the functional form by super- 
) (2) 

scripts and write explicitly Wir» Y x, for our trial 


) 
functions, and gh? ( 


(1, &\) for the f obtained from (4.10) by 
(s) ()) | 

using Vex, first and wy, second in (4.5). Tatgesett oss 

by parts, analogous to (2.50), reveal the following symmetry 


properties which are not at once apparent from (4.10): 


any J (a, wu) = FEC) 


(4.12) pv) ( - K,, K,) = { a { - Ky %2] : 


(4.11) implies that f is invariant under interchange of the 
functional forms of em? and Vie, » as will be verified 
below in the case of the Born approximations. (4.12) 
guarantees that the scattering kernel At ka, Ky) is symmetrical 


when computed by (4.10). 


Connection with Born Approximations. 
The simplest admissible trial functions are 





38 Ref. 18, Ode (2 S2)K 
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Substituting these in (4.10) we obtain 


(4.14) Rae Ve = Uy {(%, %) 


which is identical with the result of the first Born 
approximatione 

We obtain the second Born approximation if we use one 
first Born approximation trial function and one plane wave; 


1.e. either 


Uy -akel 
See - - <2 . 
(4.15) YK ls 
(2) “Bs axle er 
4h =u, ze o_o se Cake!) 
i |r | 
or 
(2) 
Y- i Y~ K 
(4.16) 


(1) 
“oat, = Sy 
in agreement with (4.11). Higher Born avproximations are 
obtained ina eer WAaYe 
Thus, apart from re-deducing the Born approximations 
from a variational principle, nothing new is gained by using these 


iterated functions in (4.10). 


Application of the Rayleigh-Ritz Method. 
We can use the Rayleigh-Ritz Method in conjunction 


with (4.10) to determine the scattering amplitude { (2, &)} 
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This is analogous to method 3 ( section I.2.) for calculating 


the tangent of the phase shift. For simplicity, we use the 


same aes form, { om YK, and eK, and BEES 


+ B ~-a kV 

-K, = 1 © + C, Ue aa. <a ce, u : 

fa, 17) & i i: eG 
ak,.¥ 


a hae 
where the trial. functions 


7 
Ui have the asymptotic form 


ak A 
ans) U(r) —> 42 (x x!) © — 
Kia J" 


Clearly for use with trial functions of the form (4.17) instead 


of (4.4), (4.10) must be replaced by the equivalent homogeneous 


equation 


(4.19) Try, 4 UT e2 C : € [ ua, ,) = Liv cy f( KR) 


aAa=L 
where tims is defined by (4.5). 


When (4.19) is written out in terms of the trial functions 


it assumes the form 


h 
« ° e ve ue 
(4.20) Z Pay , ge {O4, Ki) ¢, 
st 
y= 
where we may take As) = BS. « The stationary property of f 


relative to variations of the Cc; then leads to the equations 


x 


(4,21) 7 pS) = TH Lr, Kk): e, 


h 
2. As; ee a C2 wv 
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They are compatible only if the determinant of the coefficients 


vanishes, that is if 


By-Mi RL) Aye eee Ave 


Br, Ary ---- Alu 
: 2 fe : - | 0 
(4.22) A(K,,4,) 2 | ) | 
| 
Any Aina RO a. air, ne 


i( Kes ) is uniquely defined by (4.22). The fact that it is, 
in general, complex is reflected by the collision determinant 
N (12, ay whose matrix 1s symmetrical but not hermi tian. 

On varying the angle between ~, and <4, one can obtain 
stationary expressions for the scattering amplitudes, and 


hence for the differential cross section, in all directions. 


I-5. Non-Central Forcese 


It is welll Known that if We per talemeaumve een meee 
of a potential which depends only on their distance of 
separation, the relative angular momentum is conserved. fThis 
fact is the basis of the phase method in collision problems, 
whose variational formulation has been discussed in sections 
I.2. and I.5- 

However, the hypothesis that nuclear forces are purely 
central is untenable in view of the quadrupole moment of the 
deuteron, and has been replaced by the assumption that they 


are a mixture of central and non-central forces. Rarita and 
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Peimiiver have worked out the combinatorial aspects 
introduced into the 2-body problem by the existence of such 
non-central, so-called tensor forces. In particular they 
have shown that the angular momentum is no longer conserved, 
and that the eigenstates of the Hamiltonian, instead of being 
S, P, Dyees- States are now, in general, mixtures of such 
states. Most important, in the neutron-proton system at low 
energies, is the S + D mixture. 

In our method of section I.4., the conservation of 
angular momentum did not have to be assumed and thus it 
applies without change to collisions involving tensor forces. 
This aed deals with all angular momenta at once. It may 
however be desirable to treat each eigenstate separately. 

This can be done an a way related both to the three dimensional 
procedure of section I.4. and to either method 2 or 5 of 
section I.2. 

As mentioned abové, it has been shown from the trans- 
formation properties of the Hamiltonian that there are solutions 
of the Schroedinger equation consisting of mixtures of S and D 
states. These solutions are proper functions of the Hamiltonian, 
H; the total angular momentum J; its z-component J, 3; the total 
spin S$; and the parity operator, P.j In the scattering of 
neutrons by protons the incident plane wave includes both S 
and D waves in which these operators have the same proper 
values. The amplitudes, Bs 4 Oy of the outgoing S and D 
waves are related to the amplitudes O<¢ , Oy of the incoming 


waves by a unitary, symmetrical scattering matrix S. Thus 
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S+ of 


ol 5 
t : " ( 
fo. 2) Se on 

° " on $2.7 ° 


We know of and “). To find the differential cross section 


\\ 


(5.1) (> 


where 


R 
(t 


we must determine the elements S44 which give (3s and ‘S 

by (5.1). Actually it is sufficient to determine two of the 

( complex ) elements S-- since due to its symmetry and unitary 
character, the matrix S involves only three independent 
constants. | 


Consider the expression 


(5.3) Tie \ 4" [aran-V) 4? (ety) 


c'\ (v) 
where v and 4 are the solutions of the Schroedinger 


equation, 
—v ra a O 
(5.4) (Wa KY i 


which have the asymptotic form 





39. See ref. (18). 
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(2 emt (-) 2? 
yom {a oe age fe 


a Ls (2) AW 
+4 8 + (35 — = {Op 


eS). Sy) (t) ~AKA (\ 
Y —™* Na, — =y Ps — — | ag 


(') 


(t) —a KN aoe 
{ Ao — a (> * = oF 


here Q, 3 QO» describe the normalized angular and spin 
functions of the S and D wave respectively. Now let us use, 


in (5.3), trial functions whose asymptotic variation is 


(v3 AK : (Ee 5 aaa 
Sy —_—? ds aa € 8 (5 — {FS 


(5.6) 
| ut) cu) ety © 
a 365. — metth + §Pp a 2 
Then 
| a PCa ~) mo «4 
§Ji, = Vy = oY ~ Sy" 2 Vas 
(57) 


(2) (i) 


= Al a lest e as 56> J 


Hence, to the first order, we have 


(1) (2) 
| ey : 2 as (Ss, tine a > b>, An. ek 
(5.8) | | 
| (rv YW), 
- & | ce: (>s + ado |. 
A 
or, in the matrix notation of (5.2) 


“uM ) (2) LK ey 
+ ae ae e ae 


A 


ze M.S . 


e 
A 


(5.9) 








: << {a = 

1S ~ his 
ea fy Ee | 
7 J) = - i= % 
eee = 
De 
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This enables us to obtain all matrix elements S. e For 


(i) 
example, to find Doe we use trial functions with &@p =O , 
cul 
$ 


Su S$ 0 (vJ 
Tae (V “= ») _ ete “. Sy 
sto) (as J(s) Sa 7% Ppa te 


and an application of the Rayleigh-Ritz method will yield the 


=D in which case the right hand side of (5.9) becomes 


required matrix element. 

Let us note here that in this way we obtain a stationary 
result for the differential as well as for the total cross 
Pection ( The stationary character of the differential cross 


section is lost in the method of proper phases, described below.) 


The Structure of the Collision Determinant. 
Suppose wé want to calculate the element 31 by the 


Rayleigh-Ritz method. We use as trial functions 


(J (1) AKA AVA 
es 2 g WI yt 
V4 =H a a +05 m a s, i 45, UL +... t Ls 


AWA a) y 
+4 5 tat a 4a, ren ae -} Oy 


, (v]} ~AllA 
(5.11) vA =| le + +s, Knute — +S ae i ae mee) ‘ 


a 


Ot M4 Ov) 
+4 a jak oe WY, mie Wy, +.- (Ly. 





40. For calculating Su we need not use different functional 
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where the radial functions nee and a. vanish rapidly as 
nN—»> co. Thus each trial funct#on comsists of tam > and 
ee part.” Using aus func ti ons (5.11) ein fe. 9) we find that 
the tensor character of the potential introduces elements 
Bonnecting all Coefficients oF Ye with those of Ur ; 
The unknown §, 1s introduced by the right hand side of (5.9) 
and occurs only as factor of et Bm | 

On using the stationary property of Si » we are led 
to a system of equations precisely as in section I.2. These 
are compatible only if the determinant of the coefficients, 
hy « » vanishes. In the present ofee Ay has clearly 


the following, symmetrical form: 





S D $ D 
Nt aN ee 
Ws) w) 
Ks Xs 
l gees A 
le. a yh =X »Y a X 
| | : 
x “mix on (\) 
© oe) (0 | a Same a Y, 
° x us x u 
| | : D 
ly aim x 
by = who ie agen 
ty X44 Sy xX x XH 
| S 
Xx x XX x | (2) 
0 Ys 
x x % * Dd 
ras 2 x « | | ek 





a) (2) 
forms for Ye and Wr . We use them here because for 
calculating off-diagonal elements of S, different functional 
forms obviously are required. 
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Because of its symmetry bn (Ou) » where D. is the 
large upper rignt hand block of Dy . Su is therefore 
determined by the linear equation ie =OF Ba. comprises 
two blocks of elements connecting S and D functions. These 
describe the interaction between the S and D state. In the 
fim’ting case of vanishing tensor Toreceswthe, Vecciemzer] 


i 


\ 
and ie assumes the degenerate form 


ie D 

nn, tana TN, 

x4 954 % OO o 
S 

i 

ris) OW soles x 0 0 
0 o * w% D 

0 Oo <2 a 


Hence Ss is determined from the equation 


S 


AyD @ 
(5.14) S 


cl 
9 


2 2c 


whose elements no longer involve any D functions. This 
expresses the fact that the S wave is now decoupled from the 
D wavee In fact (5.14) has the same structure as (2.14), 


which was arrived at by considering only the S waves. 










a re 


rr my 
- a 
oo a) arse a 


ae 
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Proper Values of the Scattering Matrix. 
(x) xX) 
The unitary matrix S has two eigenstates s ( and 
with the property that | 


(5.15) | Be = Ste 2 ia) a ( wt 6 ti neah ) 


A= lle 
(.) 
The x are the proper Values of S and, because of its unitary 
character, may be written as 
iT) Dany 
(5.16) x  ~-€ 


—_ 


: XO. - 45 


Thus the eigenstates have the asymptotic form 


(1): (1) grAka Zing Gam 
oe ~ €& ’ () 
+) ~ AN i x e 
(5617) e a! 2 = 2 git 
, D Zz ) tho 


_ 1: ole (Kr +4s)(a'S Ns + a's Op) 


: . 
(5.18) Qu — > -hie 18 : sin(ine yn) (as O54 (5) 
) | 
q 


and for that reason they are simpler to deal with than other 


i | 
and pare thus standing waves with real radial functions 


states which are, in general, complex linear combinations of 


\ 
qr and cpl) « Furthermore semmneen™ has shown that the 


41. J. Schwinger, unpublished lectures. 
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-total cross section of the S + D mixture of a plane wave can 
be expressed in terms of the two constants fe “Wx only, 


namely 


(5.19) 6 <= —, = {ae “i 4 Sing ) 
It is therefore of interest that venW x and tan r can be 
directly obtained by our method. In fact this problem is very 
much akin to that of finding ‘tan " in the case of central 
forces ( section I.2. ). 

We could proceed from the most general equation (5.9) 
but the following approach is slightly simpler: 


(T) a) 
The correct x ( and. both have the form 


Sim UNA 
= fa, A e332 TEA 4 agus ti... f D 


+4 2 Stay KA + £. Ce KA a 1G an 

‘ A 2. i + 3 a 4 () 

and wy ma be taken as real ‘tnd % , 2 ’ 
zs ui ) a, rman 5 
Consider now the expression 


feel) ~ = a (v4 2 v) + (ot v) 


} 
where we admit only such iia Scie » whose variations have 





(5.20) 





where the aes bys Uy 


the asymptotic form _ 


oY Fa, a + Sa, om 4 fhe 
" .d N 





A 
(5.22) 


+h 08, Conn aa Sh, POY Oy. 



















a4, ot 
—=_t in. 0 § @ how =e 
7  S P= 2 => © Gh (eam 





_ OS eed 
7 — eS = @ ie el 
'® —a~ © oie 

jis ® : ,e =>» >, 
SS —= 


See 44 ss. © 




















nn 
——, 

























Te es . 
a = = 7 

| ao © 6 
OE ee eee . 
- rh al yle - 
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Proceeding in the usual way, one finds 


“st a 7 ea. 5% 


ial or wn 5(2)} 


Therefore, corresponding to (2.23), we find 


4} 


(5.23) 


\\ 


fect) a Tt w( A, a, + ££, ) = Kar +2) Xan 


as a stationary expression for tan Load - 


When we substitute a trial function of the form (5.20) 


in (5.24) we obtain an equation of the type 


— : vv 
oe eo) >. Vas Q ae + Qs ; Ay + Rin :4;| = Q a2) tan, 
+) J . ' 
. Proceeding in the usual way to a set of linear equations and 


their compatibility condition we find the following 


determinantal equation for tan 


S D> 
Tn ee 
fy - Klan (freee Qi, Os ae 
. | 
Fa Massey Q., Q,,..0.. S 
Tt. ; 
ey D2 bh ed bt a ele 
Qy Q,, eae Ryo Klan. Rie eee 
Qe Ce .. IR, Ki d 


@eoam ea 


— et : 
4 
% 


1 \ 





a mes 


a a 
7 — yy 
fe — 


























erin ‘oe ah 
gy es he 





— _ 
~ i ws rei 
ae _ aa | 
Peet © - ames 


> —- © 












- . 
See! a 





its 
' 










ae 





| - | ee ee ee 
_ FE OT eee 


——_ —_—— 
i al ea 





: 
| 
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Clearly this is a quadratic équation for vere » yielding 
approximations for tanny and tan Nz ° The elements Qi5 
arise from the interaction of S and D funetions. in the absence 
of tensor forces they vanish so that A factors and we are 

led back to two independent problems for the S and D wave 

( compare (2.14) ). 

From the two solutions for tann we can find the total 
cross section by (5.18), whose error will be of the second 
order if our trial wave function has a first order error. We 
can also find the differential cross section to the same order 
of accuracy as our trial function by solving the linear equations 
derived from (5.25) for the QQ: and L; » thus obtaining the 


(1/ a). 
and and then taking that linear 


wave functions ’ 
combination whose incident wave agrees with the incident S and 
D waves contained in a plane wavee To obtain the differential 
cross section to the second order of accuracy, we must fall 
back on the more general but more complicated equation (5.9). 
Numerical Illustration. 

By using the method of proper values of the scattering 
matrix, we have calculated tan “ry and tan 3 for neutron- 
proton 7S 2 > D, scattering at 2.82 MEV. We have used 
the constants of Rarita and Schwinger : to have: a comparison 
with Schwinger's numerical integration “ for the same energy. 


Thus the following interaction potential was useds 


| ain 
(5.27) Vo= a4 4(+) + 2. } tr) S12 





OO 


where 
Re 2-68 See 
f(r) en ra <a 
(5.28) = O nr va 
g(r) = ~215 203" 5 Sew 
= O 


A> 


\ 


and oT. is the spin dependent tensor operator, 


(5.29) SG = 5 [s1- © ) C4) 


NY 





PS il 


zu 


| a 
Mm is measured in units of AR = wWecO x Bo Cm. 
Schwinger’ - has shown that the Schroedinger equation 


leads, in the present case, to two coupled equations for the 
radial functions u (1) and n[a) of the S and D wave, viz. 


” 


(5.30) a | ui+ D { {« +9) = 0 


. 
ieol) | — 


an? ae ae “*) ere A (gu hw = 0 . 


\ 
where x and 9 are defined in (5.28) and le aie 


0 
is simpler to base our variational procedure on these 


equations, rather than to go back to the original three 





- > - 
_ i 
— ; = 


s 1h, F i 


;2 7t 





a *.- - - 
o = © 
= a 


” -_> - 

dl > = 

=> a 

= —-?> 's 

‘ 
‘ - 
—-_ | & => n 7 
= pia - 
—_ ie —_ «= -— «= —_—oo oe 
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dimensional Schroedinger equation. We define 
rama Av | 
= je |[S4 wey u +t O(fo4yu)| 
ol” are 
Tw | (= tee Ny ¥(geths)| } ai 


{tf we use trial functions with the asymptotic form 


(5.32) 


A ——> As( <n Un a tor vn | 


5 O35) 
we find 
= — Sy E = ) a _ a 
gf Cu = A ae +[4 5 des ur Fw) 
UDO ) 


2K ( Assam) ah Fd tony”). 


Now the exact solutions which we are looking for have tan Ns = 


= tan 0 = van le Henee the equation 


| D4 x (A? ten Ma + Ay Aa ue) 
(5.35) | : 


LAS 4 Ap ) Hany, 


wheréd "t" indicates the trial function as usual, gives a 
stationary expression for rene ’ 


We have used as trial functions 







, _ Ss i - 


' —_ on 
Pare My lil 
7 ow Pye oe 





h 
& 
Uy, = Cr, nr 
. 2 & h=l 23 
(5. 36) «ft Us | 
k+e 


I$) (d) : 
Ay ’ Ay ’ tan, and tan Nie were determined from the 
requirement that at A. = 1 the functions (5.36) pass smoothly 


mato the free particle functions 


UA = Ag | sim Ka tawny so | 


(5.37) > | 


. Ay (- jr(er) + Rome nr () 


As checks on our calculation we have at each stage computed, 

by the variational method, the S and D phase shifts when tensor 
forces are absent and compared the results with the exact 
values which are easily obtained for the uncoupled states. 


The results of our calculations are collected in Table II. 
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Table II. 


S + D Scattering with Tensor Forces. 





Uncoupled Coupled 





0.3070 









Correct Functions 9223514 0.35069 








It 1s seen that the last pair of trial functions gives 


good agreement with Schwinger's numerical integration. 


Ie6- Neutron Deuteron scatteringe 


So far we have been concerned only with two particle 
scattering processesSe We shall show that the methods which 
we have developed can be extended to many particle collisionse 


In this section the simplest collision of this type, namely that 
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between a neutron and deuteron, will be considered. Our 
procedure will be modeled after me thod © in sectvren Ieee, 
but method 2 is equally well suited for many particle 
collisions. 

The wave equation for the system of the two neutrons, 


ot" and "2", and the protons oe) omer, (0.2) and (0.5) 


| (601) I T+ V('z) + u(23) 4 Vat) / 4 =(E, +E) ¥ 


where T can be written as 


as [2%,'+3 Vy) . 
petal \ 2 \ 


(0.5) 


To remove the cumbersome factor in T, we make the following 


definitions: 


\ 


D ir 
yf. EM | va) 
Ke OW EL 





(6.2) 





42. For the definition of the relative coordinates, see (0.3) 
and (0.4). | 





58 


Thus 


63) De ~ (U4 5%") 


We now define 


(6.4) iv =¢ = ee U 


- 
NU +t En 
so that (6.1) can be written in the form 


Zt 


Now consider the expression 


| x 
(6.6) LT, = ( ae (ARE (An) (An) 
Vi 


where V is the volume bounded by the five dimensional 


a 
surfaces, 


a (0) 
A 3 A, = & 
(6.7) 


= (r) 


tv 
> 
Y 

Ww 
A 
D 


Clearly 


(6.8) -Sle= \ { ¢{-p) pF TF (-0) ¥} ry tg 
V5 | 


We aim to convert this into a surface integral over the surfaces 





09 


e—('} -~(t) 
2. and Digs and therefore express D symmetrically in 


terms of i and ¥, , namely 


2 Sox VV 17 
.* D —— (Vr — Vy, Vy, + Vr") 


&4S can be easily checked. Hence §L is the sum of three 


morms , 


(610) $L = Rare 


where 
R=} (2° woe -5F vo F*) ny en 
Vo 
feat) Lz 7 Ae, Vy, Si IL Y;. Vr & Ey %,)(ar.) 


\t 


[E* Ty Fe -d¥ Vr Y” [aty, Y(drs) 
Viel / 

According to the exclusion principle Yb must be 
antisymmetrical in the indices 1 and 2. Hence if we insist that 
our trial functions have the same property, then 8 will also 
be antisymmetrical so that M = Ke K and L can now be expressed 


as surface integrals by Green's theorems 
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Ke) W(E*H, SE -F¥ TE) Un)en) 


(6.12) Vi 
ot a a ae \ A 
Us P (Y Vy, O4 - 84 U4 de Js, (lr, | 
eG 


: W 
where 7, = ry /A, > i a 43 the three dimensional space A, €<Q 
and sy is the two dimensional surface A,=Q. Thus, 


the integral (6.12) extends over a portion of the-five 


= (t) 


dimensional surface 2 4 Similarly 


b=) su, (¥% v S£-38F%, E*)dy) ft) 
Ys 
(6.13) cy | = Wy, (¥ i Vy, fy-3¥ Vi, ©” J lary) drs) 


Vo 


al 


es ~~ A , 
LE * VIP - FE Vy E*) & ds, (de) 
0a Sa oa: 
since, because of the antisymmetry of VY » the two parts 
of L are equal. 
We now insist that a is so large that the following 
two conditions are satisfied: At ry >A , the wave function 
has split into a function (9, )deseribing the deuteron, 2 and 
5, and a function F (4, )aeseriving its motion relative to the 


neutron};i.e. 


WV 
a 


ay Y= LS) F(y,) +e e 





Gal: 


where ¥ represents a negligible contribution from decaying 


modese And further a must be so large that 


WS 


Q. 


(6.15) cp | $\) = negligible, Pi > 


We restrict our trial functions to those which satisfy 


corresponding conditions, viz. 


(6.141) rs = Cf+ (9) reir) + Ve y AZA 
and 
(6015") 4 (9,) > negligible, F = = Q. 


We assume that both <p(9,) and Ps (9,) are normalized. 
rE (") will have different forms depending on which phase we 
are dealing with. We shall assume that its angular part is 


normalized to UF. ‘Thus for S - scattering, 





Su, WA 
F(r,) = 7 5 + ®& hes. 
| A 
(6.16) } 


B fp = dan y 


and similarly for higher angular momenta. 
If we neglect ¥ and Xx » the integrands occurring 
on the right hand sides of (6.12) and (6.13) are expressed 


most simply as functions of me and Y ( see (6.14) and 


————_— FT Ol 
i = =& © Yt 
— : =i (ei — —_--< 

—_ op a 2 & 
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(6.14') ) so that it is of advantage to use these coordinates 
as varlables of integration. Therefore we transform from 


Vy, » ¥, to y » 9, according to the equations 
(6.17) 


which are a consequence of the definitions (0.3) and (0.4). 


It follows that 


(6.18) (Ar, \ (og, ) >] “(ary (Are ) 


Vn A i 3 Ve. Jy, a Vy, ) ,, 


(6.19) 


Tn) Ve ee 5 Vd, 
Now for Ay= Q > Y has the form (6.14) and & is 
(6.20) §%= g(r) 5 F (rn) -y dy ( ?) F(r }. 


Hence the gradient terms in (6.12) are 


(6.21) Ve = z cp (¢,) F [n) ¢, + ¥ (¢,) F(a) ¢, 





Vr 6 = < g'L9,) § F(n) ?, ¥ (9,) SF (x) r 


(6.22) 


+5 Sy'l9y) F(a) poe Slo) FCG) h 


A a £ j \ “i ie Cr } 
where Y, = an [> 9 Py —s fi As. and (vr, J t 
denote the derivative of F(vr,) and JF(v,) along the outward 
normal. Let us first integrate over the § \7space, keeping 
Y, fixed. The integration extends over all points with the 
fixed ¥, and fL{4Q ~- Since N,=Q , it follows from 
(6.17) that at all points outside this region 9, 2s OQ: oe 
hence, in virtue of (6.15) and (6.15!), we may extend the 
integral over the entire 9, -space. Furthermore, since ce (P.) 
and S¢1¢\) are independent of the direction of ei we can 
first carry out the integration over the angle between Y, 

A A 

and F » in which all terms involving yY, -f, give zero, 


-& A a3 
while those involving ,-V,; give Ut. Hence 


aes a 
= (=) | as | ui ey lel. 


St el Le 
(6.23) \ ¢*(¢,) Eb? (1) §€'[n) + 4(9) 5 yl¢) e* (nF (a) 
_ gt (9) F™ Un) SF(n) ~ gts) Sele) Fl) F*(n) | 


-)’\ ts, } e*(a) SF'(a) - 5 Fs) aan) 


\") 
Sq 











te ’ ite 
/ Te are Tie tery 
, 


— 











oe a ee . 


inf’, fy pcre hb bok 
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due to the normalization of (9) Notice that $ @ (9) 
does not appear on the surface. 
Let us now specialize to the case of S - scattering 


and assume that for AN,7A ’ Fy (rn) has the form 





(6.24) F,(v¥,) = 


where B is chosen to give a smooth logarithmic derivative 


across. vAjca-+ F (v,) was the form (6.16). Hence 





§f( a). dB. “eee 


(6.25) OO 
sri(a)=-dB x “<= 
3s0.)6—léuthat 
4 a 
, 3 ~ >) ied 
(6.26) K=-|2\ Rl te 


: — 
Turning now to L we observe, by (6.19), that Vy, acts 


only on functions of lo,| and that 
- A 
(627) Vi ° { Punelion 4 (9, ce “a : 


a A 
Hence the integrand contains o. Vv, as factor and the 
A 


A 
integral over the angle between 9, and ri vanishese 


Therefore, by (6.10) and (6.26) 
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| 21 
(6.28) | oS. - aK AGB 


so. that : 
c 
> TKAR 


(A 


oe ‘1 
tba + TRA BS 


(6.29) 


= a R N Hany 


is a stationary expression for tan ° 

In the derivation of (8.29) we have me t some difficulties 
which did oe Ns in the two particle problem: The neglect of 
X and y in (6. 14) and (6.14!) and the assumptions (6.15) 
and (6615") introduce inaccuracies for any finite value of a, 
even if the interaction potential has a sharp cut-off ( square 
well ); this is in contrast to two particle collision, where 
the wave function has accurately the asymptotic form for all 
distances larger than the cut-off distance of the potential, 
and where assumptions of the type (6.15) , (6.15') do not 
enter. 

This difficulty can be avoided by setting qe=oco but 
then another difficulty is encountered. We have in our 
treatment neglected all second order terms. However, unless 


Py has the asymptotic form 


(6.30) Vy imac Fx (1) 02%, > —p 


‘ _¢(5 92) Fy (n) as, > 
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where (3) is the exact deuteron function and ts ive 
exactly sinusoidal with wave number K , the second order 
terms are appreciable over an infinite portion of the 
configuration space, and their ite cates far from being 
negligible, makes I _, infinite. Thus suppose that our trial 


' function tends to 


(6.31) ty ad | lf) «dele¢ Fy ln) » 2 DA 
instead On (6480), where Fy is sinusoidal. For large No 
fe.52) Fr => (Fetau") a ee Goud 

where H -p is the deuteron Hamiltonian, 


f6.53) Lt -p E- -2 V4 V23 


Hence, from the region A, 7 2 , say, in which (6.31) and 


(6.32) hold very nearly, Ip gets the contribution 


ae | {lel +590 | Fe (4) 2 


_ Gur ey) elad+do(%) | J Co) @ed en) 


—*§ 


e 3 
s {= ) (Ep ~ Ey ) 


IF (ry dj ¢ rn). 


ok 


2 





67 


where 2 is the expectation value of Hp in the trial 
state Py ae 3 Pp - Clearly Qj although a second order term 
( since €p is stationary ) is infinite. 

We are thus faced with the following dilemma: 

(1) «If wuse only an approximation for ge we must 
take a rather small a or else one effect (6.34) which is 
infinite for infinite a becomes important and our neglect 
of second order terms is in eee However, the finite a 
introduces an inaccuracy. 

(2) If we take a carries we must use the exact deuteron 
function in the asymptotic form of 4 which may be inconvenient 
for practical purposes. ( Such a precedure would bé nearly 
impossible for more complicated collisions. ) 

The way out of this dilemma is suggested by (6.54), which 
shows just what the infinite contributions of the second order 
terms are and hence permits us to subtract theme. This is done 


by replacing 1a, (6 +66) fu Ale defined by 


mit 


cl 


( Ep [Avi )(d rp} 


Vs 


where 


V+ Ep - €p 


9! 


«\ 


(6.36) 


_, ae bo r | L- 
= Vy, = 3 V9 a Ui, ~ U2, ~Vs3, tk +E p 
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a . 
Thus To depends in two ways on the trial function: 
(1) Through P whose form is determined by the asymptotic 
deuteron function. 
(2) Through the trial - itself, when substituted in 
(6.35). 


We shall show that 
l 27 oO 
e T te oe ee 
(6.37) wis 7 TK ADS = IKAS +h, 


v 

where (\ is the integral of second order terms which decay 
exponentially in all but a finite region of the configuration 
Spacee 


Let us denote 


) ge Pe (dv An) 


cs 


Tx (4) 


AX 


6.38) l 
( \ o* Pg [ dr) (An) 
Voy 


Then, vy (6.36), 


T,' (4 


Ta (Ya) = Ta (Ga) + G!-eo)f 1 Selly en) 


- = Ta(¥) 4 STa(¥) +T0(S¥) 
sf (e'- Ey) | ZF ae.) (dn) 


Now Oo 


(6.40) Ta(¥) -0 









eet Ai ie “F 
; a fier | 


; pehij ae ;@ 
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and, as was shown above, without neglecting second order 


terms we have 


(6.41) sr, (Vv) —— rr A(B,-8) + fa 


where lim = 0 
Q = 
Soe ) 
To estimate the last two terms of (6.39), we decompose 


Va into 


(i) 
es AyAn S$ AKA 


yi) 
(6.42) a: A AANA, PR 


(Y A 
where b is so large that in V ae YY and Bt ee have their 


asymptotic forms | 
(6.43) Y = alg.) E(x] : Y, = C(¢ (F.) +5 (r,) 


or the corresponding expressions with 1 replaced by @. 


Clearly the last two terms in (6.39) are both of the second 


C 
order, so that their integral over the limited vo is of the 


(2) 


qa we shall find that the- 


type Re. On the other hand in V 
last two terms, give zero owing to the integration over §, 


) 
( or G2)> Thus in the region of ia where Jt, >t, they give 





Noe roe . (ep! -ey ) (Fal Y(Ag,) 
J Gr) P(E y-¥) + (eones ) al} apy 
(6.44) ANE Fy"~ @ F¥) (V5 + BS Sy - Hy) (pe F ~¢F) 


| (Eo ~€y) Px” Fel’ Site) 
“FPS ge (yves) oct 
= 0 


Hence 


(6.45) Es(s8} + Ver =e) Nz(" = No 


From (6.39), (6.40), (6.41), (6.45) we now find 


t.(€2) =~ we (6,~8) +)* 


which is the statement (6.37). 


This means ae ee calculated from 
} L 
(6.46) Tati wK as, 2S vn ATTA a 


is stationary if we admit only such trial functions whose 


asymptotic form is 


ean Ly — 4 (9s) \ erie i +0, keane 


Ni 


we 


| - a 
i” ¢ SOLE Tee Ha sit hs 
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Thus the exact form of the internal wave funetion and the 


binding energy of the colliding nuclei are not required. 


Numerical Calculations. 

We have made exploratory calculations on n-d scattering 
in preparation for a more extensive investigation. For 
simplicity we havé so far only used trial functions of group 
structure, but we remark again that this limitation is by no 
means inherent in our method. 

We have assumed an ordinary Wigner potential 

pews 


oo 
=~ Ve x 


a 


(6.48) V 


L ~3 44 3.) 
with V,=- 72.00 me and a = 0.8 in units of A,g= 2.80 x 10 cm. 


These constants give the correct binding energy of 2.18 MEV. 
for the deuteron. 


As trial functions for S-scattering we have used 


(6.49) f = Gx (9,) T., () ea t Ox ( ¢2) +, (rs) Kr J 


where we have set 


~059,~ 


Px (¢,) =~ O0-4237] e 


(6.50) | Fi. (rv) oa Nae Coan A ) ?) < Q {Wwe 0) (,2) 


) 
( cot d) 
43, These constants are almost identical with those of Motz and 
Schwinger, ref.23, who however used an exchange potential. 
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(6650 cont 'd) 


Ti = Fa a a(r) a Quote Scaltindnng 


Kh = 6 | aly) afr) La) n | 
al o(adels)-2 plval2) 4(3)] Doble seating 


42 2X | 
All distances are measured in units of No - Lx (9 is the 
best approximation of its type for the deuteron function. 

The polynomial approximation for E (r) has the wrong 
asymptotic Poni. Hence it was assumed to hold only inside the 
region AN <Q . | 

We have calculated only S-scattering at zero energy. Hence, 


asymptotically 


1 @! 

(6.50') cam eee A + — 
; ! 

where A and B are chosen so that the polynomial form of ¥, (r) 

passes smoothly into (6.50!) at n=QA. 

AS we are limiting ourselves to a finite region, the 
difference between I(a) and I'(a) is negligible ( see (6.35) 
and (6.36) ). The calculations to be described were carried 
out by using I'(a) throughout. 

We have calculated limit ( tan }/ = Xe, Sor quarter 
and doublet scattering, chien mite cam the total cross section 


in the limit of zero energy by the equation 





13 
Lu ev uv | 
(6653) § = HU [ S Ka a = Ka) 


where q and d refer to quartet and doublet scattering 
respectively. 
In the case of doublet scattering, the convergence 
with increasing values of W was slow. Since in the case of 
two nana Scaue eine method 2 gave somewhat better results 
| than method 3, ( see Table I, pe. 30 ) we have calculated X y 


not only by (6.29) but also by the equation 


wy Tee tt wot (LF La) — Ela) F'(a)) <0 


which bears the same relation to (6.29) as does (2.32) to 
(2.23). X is determined from the logarithmic derivative L 
according to the equation 

a’ L 
R252) ke ee oe 

al + 
However, the convergence was only slightly improved. 

The approximate values of X, obtained by our variational 
methods depend both on a and n. For large enough n we expect 
the following variation with a; | 

(Al) When a is too small the expression (6.29) for tan 
is incorrect. Hence X should depend strongly on ae 

(A2) When a is sufficiently large for (6.29) to hold, 
X should be fairly independent of a until 
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(A3) a becomes so large that F, can no longer 
approximate the true function over the entire interval. Our 
variational principle than breaks down because second order 
terms become important. Here again we expect a considerable 
variation of X with a. 

Thus we expect to find an extremum in the region (A2) 
and the value of X at that extremum may be taken as our 
approximation for X. 

As we increase n we one 

(BL) Gradual extension of the "plateau" (A2) towards 
iarger values of a. 

(B2) Increasingly smaller differences X ( a3 n +1) 
- X ( as n ) for large enough, fixed, ae | 

In an approximation method n should be increased until 
(Bl) and (B2) are so well satisfied that further increase of 
nis not likely to change X appreciably. 

In virtue of the cholce of our potential and trial 
function, all tec rete occurring in Tq could be analytically 
evaluated. Our results for X q and X gh given in Tables III, 
IVa and IVb and plotted in Graph 2. fThe best 3 = parameter 
wave functions i (r) determined in the usual way from our 


variational principles are shown in Graph 3. 





GRAPH 2. NEUTRON BEUTERON SCATTERING 
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Table III. 


Quartet Scattering : X 4 (asn) /Ao- 


( Method 3, (6.29) ) 





Table IV a. 
Doublet Scattering : Xy4 (azn) /% . 


( Method 3, (6.29) ) 
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Table IV be 


Doublet Scattering : Xa (a3n) / Ao. 


( Method 2, (6.51), (6.52) ) 





Interpretation of the numerical results: 

A. The convergence of ae much more rapid than that 
of X gs This is due to the fact that F(r ) is more irregular 
for doublet scattering ( see Graph 3 ) and therefore requires 
more terms for adequate approximation. 

Be The variation of Xq with a and nis as expected 
( pp. 73 and 74 )3; X , seems to behave in a similar manner 
although the number of parameters is not large enough to 
exhibit a proper lene. 

C. To determine the total cross section we must select 
the best values of -X q and X qd from our tables. In the case 
of X 4 ( Table III ) we have evidently come very close to the 
limit which can be reached with group structure form. We see 


that our results depend markedly on a for values smaller than 3. 





er 


This inéigates that the interaction 18 appreciable  cugaro 


about that distance. We therefore choose Xq = X (432) = “1.835 Apo. 


The choice of X y is more difficult. At a= 3 en all be 
sure that the interaction is completely negligible; on the 
other hand at a = 4 we are still. quite far from convergence. 
Let us nevertheless choose X,=z X g (432) -~ = 1.90 rn, ( see 
Table IV b ). This gives, by (6.53), S= 3.39 were. If 

we estimate the limit of X4y which could be reached by taking 
Percer Ny as Xy = = 1.50 p,we obtain 5 = 2,94 barns. These 
figures are in qualitative agreement with the experimental 
value of 3.2 barns an | 

Let us summarize the possible inadequacies of this 
tentative calculation. 

(a) Wigner forces may not be adequate. to. describe this 
process. 

(bo) Our approximation for #(e,) ( see (6.50) ) may 
not be sufficiently accurate. 

(c) The group structure approximation may not be. 
adequate. | 

(ad) The Schroedinger theory may not be capable of an 
accurate description of this process. 

Of these shortcomings (a), (b) and (c) can be remedied 
by more extensive calculation of the same kind. When this is 
done, it will be possible to decide to what extent (d) is true. 

In future calculations the following points should be ( 
given consideration: 


(1) Terms allowing for the polarization of the deuteron 
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should be introduced. 

2) Tt may be of advantage to use trial - F's of the 
correct asymptotic form. The reason for this is the following: 
If we use a polynomial expression for F (f), we must include 
a large number of terms since, due to the loose structure of 
the deuteron, the region of interaction is very large 
Bi Seece NOnus cm ). This necessitates the calculation of a 
great many matrix elements, which is probably more laborious 
than the numerical integrations which become necessary when 


functions of the correct asymptotic form are used. 


I.7- The General Case of Composite Collisions. 


The method used in the preceeding section can be 
extended to all collisions of composite nuclei in which triple 
or multiple disintegrations are energetically forbidden. The 
analytic aspects of the problem are quite analogous to those 
encountered in n-d scattering. However the combinatorial 
sice of the @eneral probibem is siomewhat complicated and no 
Aewasled description of the general case will be given here. 


We shall merely mention a few points of interest. 


Modes of Collision. 
Our picture is the following: We have a number of pairs 
of nuclei colliding and a number of pairs ( possibly different 


ones ) emergin: after the collision. The groupings of nucleons 
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into pairs are alternative (not independent) processes, 
each pair comprising the same nucleons. There is, then a 
"pepion of interaction" in the configuration space, where all 
nucleons are close together, and an "asymptotic region", 
where the nucleons have separated into pairs. 

On a hyper-surface, )_ » separating the two regions, 


the wave function has the form 


Grea) Y = 2 0. 
(i) 


where, fg, and @, ‘descuiie the normalized internal wave 


functions of two colliding nuclei and F their relative motion, 


, =< ho 
A function of the type (7.2) will be called a mode of 
collision. 
ir Y is a proper function of a complete set of operators 
(H, J, P, 43..), its modes satisfy the following 


asymptotic orthogonality conditions: 


— iG 
(7.3) \ 2,9, lf 20d ) AFY- 
y , 


44. The term "collide" will be used in the remainder of this 





section to describe both collisions and disintegrationse 
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Hort f O. and ?, describe different palrs Ce enuc lca, 
they exist on different portions of the hypersurface Ze and 
hence (7.3) is satisfied. If however they describe the same 
paiv Of nuclei, then emu @; +e} or Bi 4 es so that the 
orthogonality (7.3) follows from the integration over the 


internal coordinates of the nuclei. 


Number of Arbitrary Constants of the Scattering Matrix. 
We have remarked in section I.5. that due to its unitary 


character and symmetry, the two dimensional: scattering matrix 
involves only three independent constants. In the present 

case the scattering matrix has the same order, n, as sthere are 
aifferent modes of collision. When represented in EHS scheme 

of incoming and outgoing modes, we have proved that the number 
of independent constants is (n) (n+l). Only this many constants 


need to be determined. 


Rayleigh-Ritz Method. 
We have already applied the Rayleigh-Ritz method to two 


wo blens involving two modes of collision: n-p scattering with 
tensor forces ( section I.5. ) and n-d scattering ( section I.6.). 
In section I.5. we have discussed two different variational 
methods; One dealing with complex ( running ) wave functions 

and yielding all elements of the scattering matrix; the other 
dealing with standing waves and giving only the proper values 

of the scattering matrix. ( Only this second method was used 


in ned scattering. ) 
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These methods can be generalized as follows: 

Running waves; ] 

When dealing with a special case of running waves in 
section I.5., we introduced a wave function, Y ; 
(5.5), whose angular and spin functions were the complex 
conjugates of those of the function We - This was done 
in order to make use of the orthogonality relation, 
corresponding tO (760) 

In thee presentemore general sease, such” a pmocedurs 
is not convenient, since the wave function does not everywhere 
separate into radial and angular parts. On the other hand, 
if we take the complex conjugate of the entire romevret,, tne 
sense of the incoming and outgoing waves is reversed and the 
ck 's and 0 's, denoting their amplitudes are mixed up. To 


avoid this difficulty we use the following notation: 
(7.4) ya a Di) * 

) 
where ee and yey) have the asymptotic forms 


— (’) — iM (v] peas AKA 
se = a ¢ Gs Yi (ae _— «Op ew )s, 
4 an 7 Ay, 


(7.5) i Du De 
yee : 5 ~ ie ae KA, : (=) e <= 
4 : A, (4 A, : 


S 4 is the angular function of relative motion of the mode Q . . 


We note that by the definition (744), 


a ® a ° 
A A 3 


(7.6) a. ie (057 aie 
4] 





82 


ae 
which is the reason for introducing yr as only in this way 
can the o's and (> 's ‘remain associated with incoming 
and outgoing waves respectively. 


We now form the expression 
we) Lh = 4 ~U + 
V 


where y a ' oe the same quantum AanbeOn Vis the 
region of interaction, H the total Hamiltonian and E the 
total energy of the system. The first variation of I2, 
can be transformed into a surface integral in the usual 


way and we are led _ to 


ria oS f¥OrD FEMS EHD gr 


Ou “Orn 


i7aS ) 


> ae ve ee 
khaoe 50: ( Ki /M; 
ba > Nod. mass 
(Compare (5.7)) Therefore the following equation is 


stationary relative to variations of Y and ve . 


ia (-1) n Ce? 
(7.9) le —_—* Ke . oO - ‘ —_ ate ° ol. ° 
u\ ai e 2 = a ry Anat mis ~~ 9 K. a pl 
We would like to convert this into a stationary 
expression for the elements of the scattering matrix 


(see (5.9)). Now in the case where groups may collide 


with different kinetic energies, we must modify the 
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definition of the scattering matrix given’in (5.1). The 


amplitudes (> are linear combinations of the Q} , viz. 


(7.10) (>; = = Cc 
) 


The total flux into the region of the configuration space 
enclosed by os is 


(7.11) a2 K; | a; \"- us 2 Jers | of 5 vt 


22 Af al ETB es Nm; ‘tj 


; G 
The last line shows 4 tnat the matrix 


(7.12) S = [ es C35 ) 


a) 


1s unitary, and hence may properly be taken to be the 


scattering matrix. If in analogy with (5.2) we define 


A ig, of ' 2 I, (>) 
se ; {z, dX, ) 7 (it, Pe 
; | 





46. Ae Wintner, Spektraltheorte der Unendlichen Matrizen, 


(Leipzig, 1929) pe 34. 









: ay = ees 8 @ at 





84 


then by -(7.‘1vu) ana (7. aay, 
14) © == SA 


In this notation (7.9) can be written as 


aes pgl tsiagil"2 Oy 
r (2) a 

ee LO} a“ a 8 . = eae A 7 S i A 

: : z\ + Be A B Pi wt A 

Hence by choosing our trial functions with only a single 


non-vanishing component of A, we can obtain a stationary 


expression for any one matrix element of S at a time. 


Standing Waves: 

The modes of the proper functions of the scattering 
matrix all have the same phase shift (see section I.5.,(5.15)4t.) 
We form 

¥ z | i 
V , 


—” 


where Y 1s a trial function which consists, asymptotically, 


of standing waves. The first variation of LT now becomes 
k> 
iol?) 4 <-8) K: A. 5 Harsh ; ; 


where A: is the amplitude of sw KN: [y; in F Gene 
function of relative motion. 


Hence 
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The collision determinant, which arises on application 
of the Rayleigh-Ritz method is of the a’ denice corresponding 
to the sum of n squares multiplying tan Ae Its n solutions 
are approximations for the tangents of the n proper values of , 


the phase shift. 
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Ii. THE VARIATION ITERATION METHOD IN COMPOSITE COLLISIONS. 


Hl. Oe Pre jimi Wana Considerations. 


Scimiineen' Nas receiv, developed a variation iteration 
method for two particle collisions. His treatment leads to a 
stationary expression for vat » involving a Green's function 
and giving very good results, even with rather crude trial 
functions. In this section we shall first rederive his 
results in a slightly different way, to show its relation to 
the bilinear variational principles discussed in Part I; and 
then shall examine this method with a view of generalizing it 
to many particle problems. 


Consider first the one dimensional Schroedinger equation, 


lt 
(8.1) Re 4 i & ® £(x) \ u(x) = 0, 


where u satisfies the boundary conditions 





(8.2) u(o) =O ; u(x) > Sunwie OC 4 Vany. co» KX 


OA 0% VW CO. 


Let the Green's function G(x, x') satisfy the differential 


equation 


(8.3) s aw n° } G(x, x') sag 3 ({x-x') | 
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and certain boundary conditions to be determined later. 
Multiplying (8.1) by G and (8.3) by u, and integrating over 


x from 0 to a we find 


ct r 7 
G (x)x') ae —- ulx) re. C ( x, x') | 
(8.4) 


a 
: +a C (x, x') f(x) u(x) ole = u(x'] 


To make the integrated terms vanish we now impose the same 


boundary conditions on G as aré satisfied by u, namely 


(= [ 0, x') au 
(8.5) 
G xm xc!) —_— 2? F(x) (ciuk x + Howry. Coax) Osx > oo, 


where F (x!) is some function of x!'. The solution of (8.3) 
and (8.5) is 


CoA 
(8.6) G(x, x') - Sum KX, - Cod Ky + Tk sin KX. Stun Ke 


\ 
K 
In virtue of (8.4) and (8.5) we now obtain, on letting A 30d, 


oD 


fee7) ujx)= A | G (x, x!) He) ul») da, 


This may be regarded as a proper value problem for A, 
for a given coty , and the discussion can. be completed from 


this point of view: (8.7) is multiplied by f(x)u(x) and 
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Lategrmabed over x, giving 


F tlo[ ute] de a = a | 5 WON: ¢ () u(x) ae | : 


(8.8) 


Mf) a) Cla {Ua ls!) dx da 


where 


(8.9) Gy(x>x!) = Sink X, + CoomM Zo) 


( 
K 
(8.7) is a stationary expression for x corresponding toa 
a given cot ue cg corresponding to a given a ° 

On the other hand we may derive (8.8) as follows: 


We note that u, defined by 
o® 


(8.10) u, = A | Citas f(x ') ul') dhe 


becomes, asymptotically 


a (Sim ka =f Hom Cw ex), 
(8.11) 3 a 
dort | See x f (x') u(a') dv 
K 
0 
and thus has the correct asymptotic form and phase shift 
for any trial u; furthermore by (8.7), u,<u when u is 


the correct solution. Hence if used in our bilinear 
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variational principle (2.48),in conjunction with u,= u, 


we have simply 


tai = er Loe we RAO] 


(8.12) a (Oo (x,x') + _ Si Yer. SL. Kx' ) {Ota (oy de | che 


b 


= 5 5. ; 
) 


5I5, = O 
for arbitrary variations of u. This is of course 


equivalent to (8.8). 


Three Dimensional Formulation. 


Turning now to the solution of the equation 


(8.13) |v’ ie oN f(r) |y(ed=o. 


one can proceed in an entirely analogous manner. The solution 
of (8.13) corresponding to a wave of unit amplitude incident 


along (4, is written in the form 


| ee : 
(8.14) Wy = ete oe a} : flr vie: (v')(ar!) 


“Ge leer' 


and becomes asymptotically 


A Lt IT 


5 a7 ZW 
(8.15) Ve, —> e a + ce | ¢ =a + ar Yur iar’) 


Hence the scattering amplitude in the direction K is 
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given by 


wae (How) = Efe {4 gle le). 


Ut 


This permits us to write (8.14) in the homogeneous form 


2 Yale) = Af C(ne) PO) yale) (bY 


where 
{ aK, -7 Likyr o 4 Klr-r'| 
(8.18) G(r, r') >, =a en e <a 


Hr f(xy — 


For large values of r, (8.17) becomes, in the direction IK, , 





J | — aw or oh. 
Ee 4m 4 (kr, ,) (° ¥{ (v2) A 


(8.19) | 
Je et) vale eed 


h ; bs ti 

Hence no matter what function Yu, (r is substituted 

in the right hand side of (8.17), the left hand side always 
has the correct scattering amplitude f( 1K) )- : 
This, in fact, 1s the.reason for replacing (8.14) by (8.17). 
(t}-10 


Returning to the .variational principle, it is clear 


that if we use only such \/y which have the form 


2:20) Yul) = OW] Clue) €O') ea(ler) , 





oa 


(v) 
and Benny Le, (r) » then since {+r (42, K,) = K (Re, K,) 
it reduces to ty, Ky =O When this is written out in 
full, it serves as a variational principle for the 


determination of 4 (1,%))3 





ox Aw-V Hy i | 
i xi | Yl) fle) 2 Ar) \ Paster) fle) oh (de) 


20) — leu) {1 Gals) (Av). 


eT r'| 


t Bye (r) f [+ ) ) £0r') lr} (Av') 20 


Let us note specially that no boundary conditions had to 


“un [rr r'| 


be imposed on either Yu, oF Q ka » This is in keeping 
with the fact that all integrals occuring in (8.21) extend 
only over the region of interaction, where f(r) and f(r!) 


are not zero. 


The one and three dimensional treatments have in 
common certain features which may be expected to carry 
over into a similar treatment of composite scattering 
processes: 

Ae Both treatments involve one iteration( (8.10) 
and (8.20)) and hence give excellent results even with 
poor trial functions. In principle the accuracy can be 


incroased by itorating a number of times. 
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B. The trial functions need not satisfy any boundary 
conditions; instead, | 
a C. The boundary conditions are automatically satisfied 
by a proper choice of the Green's function. 
D. The Green's anc wen 
1. except for the Inhomogeneous term, satisfies 
the same differential equation as does the wave function 
asymptotically ( see (8.3) and (8.18) ) and 
2. has the same phase shift ( or scattering 


amplitude ) as the correct wave function. 


IIe9.- Neutron Deuteron Scattering: Six Dimensional Green's 


Functions. 


As a prototype of composite scattering processes we 
shall consider in this section the collision of a neutron 
with a deuteron. -We shall generalize the method described 
in the last section in a straightforward manner, but shall 
find that for composite scattering it is too complicated to 
be of practical use. | 


The wave equation for our system is, by (6.3) - (6.5) 
UL 
(9.1) PY — (=D -un -vrs-oy, 1H 49) G0, 


where 





93 


(9.2) D = - Wr + $ Vy, ) — (Ya 4 : Ve), 


Let us for the moment ignore the spin dependence and 
simply insist that the wave function be antisymmetrical 


in 1 and 2. Then, asymptotically the S-wave has the form: 


i ee 


a 


| Kar. 4. A 
oo — - 9) Sint ee 
L 


ee, oo Lf —? ely) 


(9.3) 


The wave equation (9.1) becomes asymptotically: 


i — J oP . Ly Veo —> (vers K™ ) +(4 ee tus 46n) \F=0. 
(9.4) 
ae v % 
Oo) —> oD ; L Fos | (Were) «(4 Ve, -~Vi3 + eo) |¥ 2 On 


According to Del., section II.8, the Green's function 


should satisfy the differential equation 
foe5) G=-4 (r-r,') J e- r,') 


in the asymptotic region of the configuration space, 


where L reduces to L, and L, for large values of fn, and Ny 


\ \ 
respectively. In the one dimensional problem L was simply 


u 
a + Ko and made no reference to the interaction 
dat * 


potential; hence a solution for G could be easily founde 


— — 
e 7 
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The fundemental difficulty in the case of composite 
collisions is that parts ar te interaction, namely the 
internal forces of the colliding nuclei, persist at infinity. 
Therefore they appear in the asymptotic operators L and L 
and hence also in L, which makes the derermiertet on of G 
practically impossible. In fact the simplest L which we can 
choose is L=P+v,,, since v,, is the only interaction 
which disappears asymptotically. Of course the practical 
solution of 
(9.6) (- D — U2. - Vay + K+ Ey ) Ce > — 5 (44- a9) d(v,-9) 
is quite impossible. 

For this reason we shall not pursue this method any 
further, but merely state that if formally carried through, 
it has the properties A, B, C and D mentioned in the last 


section and leads to stationary expressions for the elements 


(9.7) (ue? ) Ka ie \,2 ; y= 


of the scattering matrix ( 1 and j refer to the incoming 
and outgoing neutrons), \4 oe BSG to their directions). 

Another formulation which suggests itself is analogous 
to the treatment of the bound state: We write the wave 


equation as 


and introduce a Green's function to satisfy 










ee ee | 
a — ~*~ ee — - a! ¢e 
7 oe Ss a> 
ns © sad —_— 
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(9.9) (Vi, +4 Ve, ) G Fs, °3 r °,') = ~ 9 (r,-y,') 3( 9, ta 


A SOlution of this equatien 13 





i oO 1 
(9.10) G ree —_—_ oar ar L Zz " feos v 
as yaie 97% 528 0 (i (er) Sms) 
| 
where 


vi lv3-C; 


p= 193-9; | ‘ 


Since vn Ep is negative in our case, G falls off exponentially. 


(9.11) 


Hence, on multiplying (9.8) by G, (9.9) by 4 » subtracting 
and integrating over Y and 93 » the surface integral 
vanishes and we are left with 


waaay F (W951) 2) 6 (sg 5 1) (Oaronston) HoH ks 


However the exponential vanishing of G wipes out all reference 
to the boundary, from where in the past we had always derived 
the connection with the scattering properties. Thus the 
present formulation, while providing an iteration procedure 
for the collision wave function, does not lead to a stetionary 


expression for the phase shift or scattering amplitude. 


efit? a’s 


i? 
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II.~.10. Neutron Deuteron Scattering: Three Dimensional 


Green's Function. 
' Since, as we have seen,the straightforward generalization 
of the two particle procedure is of no practical use we have 
tried the following alternative approach: 
Where ¥,3 208, the wave function separates and therefore 


satisfies not only,L,Y=0 ( see (9.4) ) but also 


(10.1) ¥,—» <9 (Ver K~ ) Y (9, ) i) =U 


Hence, in partial agreement with D.2 ( section II.l. ) we 
introduce a three dimensional Green's function, which | 


satisfies 


(10.2) (Jr 4x”) C (riyr,') _— § (4-1) 


and has the correct scattering amplitude in the direction 
KG > Suppose that we consider the wave function, Wb cls 


corresponding to a plane wave of neutron 1, of unit amplitude, 
46 


incident along |4, while both neutrons 1 and 2 are scattered si. 


Then the scattering amplitude of neutron 1 along 6 3 


re wn « A Green's function satisfying (10.2) is 





46. This is in conflict with the exclusion principle but 
convenient for our present considerations. By super- 
posing the symmetrical wave function, with a negative 
sign, a solution obeying the exclusion principle can, 
of course, be obtained. 
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okt nent 
Glyn!) = Se 
(10.3) ua [n-v,'\ 
elt or 
im \ e. KY, ae Ke 1 


ae {CE 8) 
It satisfies the differential equation and asymptotic 
conditions of the wave function only where Y, Sce© , but 
not where = oO, and thus not over the entire asymptotic 
space ( see section II.l., Dae 4 
Proceeding nonetheless we write the wave equation (6.5) 


for our system in the form 


— 


-*) (Vr 4K?) YL el L¢,,%) = UL ga") Y Ru len) 


where 


= _ New , 
(10.5) - W954) = 3 V 4 Vi3 +05, +42 — Ep 


S1 


is a kind of effective potential, acting on particle 1 when 
f, is small. — 

From (10.2), (10.3) and (10.4) we can deduce in whe 
usual way that the wave function satisfies the integral 


equation 


aoe) & Ke le, r,) = Veln, “') Ul es ) F yy (0, n\(d iv ; 


as may also bo checked directly. For large rf, in the 


direction Ko and small values of 5 this becomes 





3. 


( 


ite 
t ( aK, -f ay WW 
Ky (g,5" ui (Ks et) i 4 f(x’ oo 





Ge 7 } 


+e 
eats" Ue.) Finle, nt) [ds | 


Thus again, regardless of the YX ul) (9, 01) used in the right 
hand side of (10.6), the left hand side has the correct 
scattering amplitude along - However unless 


! 
splits for large fr, into a product of the form 


(oe) af (geen) —> 4(9,) F(x) 


UY does not vanish so that (10.6) would diverge. 

A 

For in contrast to O4(r) UW is not a simple multiplicative 
factor and "cuts off" only wave functions of the type (10.8). 


Now consider the expression 


(10.9) [<= Y_ a 8 (Ve 4k ul gyn) ) F0 (dv) (de), 


Clearly, for the correct Gy! ,» £ =O. Let us now consider 
the variation of I, if we admit only trial functions of the 


form 


(10.10) Y! v2 16 (runt U ( ¢,, nv) Oy ule, x w')(de) 


=e 


where ® will be suitably restricted, and af | is 


unrestricted. Then 


avr, 


ne 


A\ 
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(Oe, 1a) eal \ \ “Kp -D) 8 Lye ~~ SY (-0) Fyfe "eke, 


As in section I.6., (6.8) = (6.13) this leads to surface 
Un 
integrals over the five dimensional surfaces 2 me | 1, ee) 


aaa ( a ir 
and a (r2=4) » The one over es vanishes because =o 


Wty 
has the correct 4( Ka a Ko , owing to its expression in terms 


\r = 
ne - ir 3 Hoy 


of a Green's function. The one over Zk 

describes only outgoing waves of neutron 2, since Yay 

contains only such waves and hence the integrand vanishes. 
| We shall show that Shey) has indeed this property 

( of outgoing neutron-2 waves ) if only Qe has it. For 

Oya and NSQ, i excecac 2/3 Be 


Hence, by our assumption, 


(10.12) Dy, os (92): {eo 


Now for 0). % 2/5 a, (10.12) is an exact solution of the wave 
equation (10.4) so that 


(10.13) hae Ro 4 Uy, | D xl = 0 


Hence 


Fv 


at 


\¢ (ry 1, Ulg, ¥ ] D yw (9, r] dy, 
-\ G (ny ni) (Sy 4 “) ® Ka (4, Nw) (Ar, ) 


qi 


(10.14) 


\oy U8) (Tr, +e) Clwy) & v) 






a Fn * ie 
ae 


pal oye ue Fel 
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(10.14 cont'd) 


| -{ \ (w4n') Vr, @ £ Dy (pn) - 
| D yi (fn) Vy, G (sn) | dS, 


But when 2 Oy a and ria & is (ep. J ‘= 0 
7 ) la ) } | 

so that the surface integral vanishes. Hence, using (10.2), 

we find 

aw, 
Db he 

10.15 y ——— : (t) D ("J — { 

( ) - 2a rye = K, els ») 


7 


Ay 


as was to be shown. 


Therefore if the trial function Oyw satisfies the 


Menai tions 
_ Ahi Jor ls.) Fla) as, > op 
(7 Qe16) D qn) ; 
| aay (92) {. <— 


the equation ) 


= | Yn (gyn) U oat ") dD ar fe ' ri) (ap \ier,) 


\) . ak iN 
Ojua:'7,) ; \Y xs (gs) le, "] C \ | ; 


uw [y,- r,'] 


e7 Wr 





Ulin!) Bam (g,4')le, dlrs) (de) 


( aor ee a) 
+ Halal mre ® Ulan) Dy 
t _ leg dienyl). 
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provides a stationary expression for {(ka, Kor 

In this way we can obtain the amplitude for elastic 
scattering in all directions. In the simple case of n = @ 
scattering, the amplitude for exchange scattering is then 
determined owing to the restrictions to which the scattering 
matrix is subject. Thus, let us consider a sub magia (Sz ) 
corresponding to scattering with a given angular momentum, Ly 
where the subscripts refer to the incoming and outgoing 
neutrons. There are two modes of disintegration so that, by 
section I.7. ( 537) involves at most three arbitrary constants. 
In fact, due to the equality of the two neutrons this number 
is further reduced to two. Thus a single diagonal element com- 
pletely defines the matrix. It is clear therefore that the - 
determination of L(g") Ke) from which all elements ae can 
be deduced is sufficient. 

Variational principles for $ (Ki) Kp} can be obtained 
in a similar way, but in addition to volume integrals, as 
occur in (10.17) they also involve surface integrals. The 
same holds for variational expressions for the proper phases 
of the scattering matrix. | 

No numerical application of (10.17) or its equivalent 
for the phase shifts has been made. The labour involved 
would be enormous. Thus the calculation for the phase 
shifts involves a six dimensional integral which can be 
reduced to five dimensions by an analytical integration. 


The remaining integral, owing to the presence of the Green's 


function must be evaluated numerically. Furthermore, since 


eee SO OS ee 
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the iteration (10.10) extends only over three dimensions, it 
is not certain how much 1t improves the trial funcr venmeone 


hence one cannot be sure that (10.17) will yield very good 


results. 
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CONCLUSION. 


The problem of ace particle nuclear col 14d Tone bs a 
difficult one. Even in classical mechanics the treatment 
of several interacting particles is in ponenet complicated 
and hence one cannot reasonably expect a royal road to the 
solution of the corresponding nuclear problems. 

As remarked earlier, however, good results have been 
obtained for many particle bound states, with a Peaaonete 
expenditure of effort and time. In Parts I and II we have 
examined the two me thods which have veen used for bound 
Swisuems, L0r eae applicability to collision processes. 

We believe that the Rayleigh-Ritz method, discussed 
in Part 1, provides a practical procedure also for the 
solution of collision problems. To be sure, the fact that 
eae trial functions must satisfy certain asymptotic conditions 
(e.g. (6.14) ) causes some complication. But the integrals 
to be evaluated in a collision calculation are of precisely 
the same kind as OCCUr in thewbandiag energy calculation of 
the corresponding bound state. This indicates that scattering 
cross sections can be obtained by an amount of numerical work 
comparable to that tnvolvean bound state calculations. 

On the other hand, the variation iteration method 
( Part II ), although yielding strikingly good results in 
two body collision problems, appears to be unsuited for the 
treatment of more complicated scattering processes. For 


we have seen that the theoretically most satisfactory procedure 





104 


( section II.2. ) is completely useless from a practical 
point of view, while other netHens: ( section II.3. ) are 
still extremely complicated and of doubtful value. Of course 
one cannot rule out the possibility that some better iteration 
procedure exists. However any integral formulation increases 
the dimensionality of the integrals to be found. Tnus, unless 
some unforeseen simplification can be made, we do not believe 
that any such method can be as useful as the Rayleigh-Ritz 

me thod. 

While this paper falls far short@irem providing 
complete solution of the problem of composite nuclear 
collisions, it is hoped that it has made some contribution 
towards paving the way for a systematic treatment of cont 


processes. 
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